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1.0 INTRODUCTION 


PREFER is a recursive orbit determination program which is used to 
refine che ephemerides produced by a batch least squares program (e.g., 
GTDS). It is intended to be used primarily with GTDS and* thus, is com- 
patible with some of the fiTDS input/output files. 

The Input to PREFER consists of five disk (or tape) files and the card 
input. The satellite Orbit File contains the nominal satellite ephemerides 
and the state transition matrix as generated by GTDS. The ephemerides on 
this file should be accurate to within approximately 50 meters (to minimize 
linearization errors). PREFER interpolates from this file at the times 
given on the Measurement Data File. It processes these measurements in a 
Kalman filter to estim^itc the corrections to the nominal trajectory. The 
filter state also includes other parameters which have an effect upon the 
orbit determination (e.g., drag, perturbing gravitational accelerations, 
thrust, measurement bia,f>es and refraction parameters, etc.). 

The measurement data types which PREFER can process include ground 
range and range rate measurements, TDRS relay measurements (range and range 
rate), GPSPAC pseudo range and pseudo delta-range measurements, NAVPAC 
range difference measurements and altimeter measurements. Processing of 
the ground measurements requires that earth motion be accurately modeled. 
The Solar/Lunar/Planetary (SLP) File supplies the data which is used in the 
true of date to 1950 coordinate transformations while a Time Coefficients 
File supplies data used for time system transformations and polar motion. 

A GPS Trajectory File supplies the ephemerides of the GPS satellites which 
are required to process the GPSPAC or NAVPAC measurements. TDRS ORBIT 
Files supply the required TDRS ephemerides. 

After running the Kalman filter forward to the end of the Measurement 
Data File, PREFER performs optimal smoothing. A file created by the Kalman 
filter is read backward in time and the smoothed estimates are obtained by 
using the recursive formulation of Rauch-Tun-Streibel . 
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The combination of a Kalman filter and a smoother should result in 
greatly improved estimates of satellite ephemerides as compared to batch 
estimation. Batch estimation is subject to errors because of errors in the 
dynamic models (e.g., gravitational). A filter/smoother which properly 
accounts for dynamic (state) noise should weight the data optimally and 
reduce the estimation errors. Smoothing will produce better estimates (in 
the middle of the data span) than just a forward filter because past _ar^ 
future data is used to estimate the state at each point in time (a filter 
uses only past data). Smoothing also tends to average out any dynamic 
modeling errors which remain. 

It is suggested that the user carefully read Section 3.1 of this guide 
and the PREFER Mathematical Description before attempting to use the 
program. An understanding of the dynamic models and correct choice of 
state noise are essential to proper functioning of the program. 

PREFER requires about 370K (base 10) bytes of storage on the IBM 360/ 
95 for a run using only ground tracking. Since the core storage is 
dynamic, runs using TORS or GPS tracking may require as much as BOOK. 
Execution time will also vary drastically depending upon the data span and 
the number of adjusted parameters. It is expected that a one day solution 
containing 2200 measurements processed in 120 sec "mini -batches" will take 
less than 2 minutes of 360/95 CPU time and 5 minutes of I/O time. Approxi- 
mately 6016 of the CPU time is spent in filtering. The I/O charge will 
almost double when the option to output the smoothed covariance is used. 
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2.0 INPUT FILES 


The following sections describe the PREFER input files mentioned in 
the introduction. 

2.1 Satellite ORBIT File (Unit 50) 

The satellite ORBIT File is generated by 6TDS on unit 21. This file 
contains the ephemerides and matrix sums used to compute the state tran- 
sition matrix ( parti als of current cartesian elements with respect to 
epoch elements) of the satellite. Although the file may contain parti als 
with respect to many paramenters, PREFER only uses the 6x6 matrix corre- 
sponding to the epoch orbital elements. 

The ORBIT file consists of three types of records: two header records 

and the data records. The file is written on a tape (1600 BPI) using a 
record format of VS, a logical record length of 6664 bytes and a bloc's, size 
of 6668 bytes. The format of the individual records is given on the follow- 
ing pages. 


3 



Satellite name In EBCDIC 

Area of satellite (km^) 

Mass of satellite (kg) 

Satellite reflectivity constant 

Drag coefficient 

imiMDD. of start date 

HHHMSS.SSSS of start date 

YYMMDD. of end date 

HHMMSS.SSSS of end date 

TIMMDD. of epoch date 

HHMMSS.SSSS of epoch date 

YYMMDD. of reference line for time 
coordinate system 

Greenwich hour angle at epoch (rad) 
Julian ephemeris date of epoch 
K-1,6 

Keplerlan elements in coordinate 
system of integration 

K-i,6 

Spherical elements in coordinate 
system of Integration 

K«l,6 

Position and velocity of satelllt 
of Integration 
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RECOnO 

RECORD 

SIZE 

•YTES 

name, description, etc. 

1 

6660 bytoo 

257-368 

OBLIMT(L) : 

L-1,J.4 

Auxiliary orbital aleoancs In 
coordinate system of Integration: 





L"l, Eccentric anomaly 
2( Period 

3| Time derivative of period 

4, Mean motion 

5, True anomaly 

6, Perifocal height 

7, Apofocal hei^St 

8( Time derivative of argument 
of perigee 

9| Time derivative of ascending 
node 

10» Velocity at apogee 

11, Velocity at perigee 

12, Latitude 

13, Longitude 

14, Height 



369-376 

0BSYW>: 

YYMHDD. of start of fitted data 
span for element set 



377-384 

OBSHMS: 

UHMHSS.SSSS of start fitted span 
for element sec 



385-392 

OBEYMD: 

YYMMDD. of end of fitted data 
span for element 


1 

1 

393-400 

OBEHMS: 

HHMMSS.SSSS of end of fitted data 
span for element sec 



401-408 

WBMS: 

Weighted KMS of fit for element sec 



409-576 

COVMAT(I) ; 

1-1,21 

Upper triangle of state covariance 
matrix 



577-584 

AZERO: 

Difference between A.l time and 
UTC at epoch 



585-592 

TZEROl: 

Time from beginning of year in 
seconds 



593-600 

DEPOCH: 

Julian dace of epoch (A.l system) 
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1 

F(fCC«iO 

Rccono 

SIZE 

■YTES 

NAME, DESCHIPTION, ETC. 

1 

1 

6660 bytes 

601-608 

SPARE 

Spare location 

ii 



609-612 

IDSAT: 

Setclllte nunber 

’1 



613-616 ! 

NBRRUN: 

Run number 

1 



617-620 . 

NEREIS: 

Element set number 

i 

% 



621-624 

150: 

Inertial coordinate system 
reference indicator. 1-1950, 2-TOD 

1 

i 



625-1024 

INDSEC 

(I.J): 

1-1,20; J-1,5 

Force model indicators 

■ ,4 



1025-1048 

mcdwf(k; : 

K-l,b 

Model Identifiers 

1 

. j 
/ 



1049-1052 

NBROBS: 

Number of observations In fitted 
data span for elements sets 

• 1 

1 



1053-1056 

NSTATE: 

Number of state partlals 




1057-1080 

KSTATE(K) : 

K-1,6 

Label numbers of state unknowns 




1081-1084 

IPART: 

Partial Indicator 
■1 partlals on data record 
■2 no partlals on data record 


1 


1085-1088 

ICENT 

Cantral body Indicator 

1 : i 



1089-1092 

IND(l) 

Orbit generator indicator 
- 1 , time-regularized ^Cot'ell 
orbit generator 
■ 2, Cowell orbit generator 

■\ 

A 

A 

■ j 
! 

i t 



1093-6660 

Spere 

• 

Ii 

•1 

2 

6660 bytes 

1-8 

DTIM 

Duany t’ime word for internal 
retrieval use (± 99999999.0) 




9-16 

ZERO 

Zero word for internal retrieval 

uae 

i- 



17-104 

GM(I) : 

I - 1,11 

Gravitational constant times the 
mass of body (km^/sec^) 

i 
1.1 ■ 

M 






nicono 


nicofio 

SIZE 


■rra 


2 


6660 bytes 


105-144 


145-384 

385-408 

409-608 

609-612 

613-616 

617-620 

621-624 

625-628 

629-632 

633-672 


673-676 


677-680 


SECM0D(36,I): 
COEPl(I.J): 
SAE(I) : 

SPARE 

IND(40) 

tND(41) 

CND(42) 

IND(43) 

IND(44) 

NTAB 

JARGl(I): 

KONFIG 

NEQ 


name, oEscnimoN, etc. 


I - 1,5 

Tins regularized exponent of 
satellite radius for Section I 

I - 1,10; J - 1,3 

IlAe conversion coefficients; 

I Indicates date Interval 

I - 1,3 

Drag coefficient tloes area of 
spacecraft surfaces 
I " 1, sphere or end of cylinder 
1*2, sides of cylinder 
I • 3, paddles 

Spare locations 

State partlals indicator 

I ■ yes, 2 ■ no 

Drag partlals indicator 
1 • yea, 2 ■ no 

Solar radiation partlals 

Indicator 

1 ■ yes, 2 ■ no 

Potential partlals indicator 
1 “ yea, 2 ■ no 

Thrust partlals Indicator 
1 - yes, 2 ■ no 

Time conversion table Index 


I - 1,10 

Julian dates which define date 
Intervals for time conversion 

S/C configuration switch 

0 *• spherical 

1 - cylindrical 

2 <■ cylindrical with paddles 


Number of variational equations 
to be integrated 
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RKORO 

micono 

SiZf 

•VTO 

NAMl, OiSCRimON, ITC. 

2 

6660 bytaa 

681-684 

NCNM 


number of to be estimated 


* 

685-688 

NSNM 


Number of to be estimated 



689-768 

NDEG(I) : 


I - 1,20 

N indices (for Cjj and 



769-848 

CiORDCl) : 


I - 1,20 

M indices (for and 



849-868 

INDY(I) ; 


I - 1,5 

I" 1, array location for drag 

partials in the Integrator 

I ■ 2, array location for drag 

radiation partials in the 
integrator 

I ■ 3, array location for drag 

potential partials in the 
integrator 





I«4i. Array location for thrust 
partials in the Integrator 
5, spare 



869-872 

IND<39) 

Transformation Partials required 

1 ■ Partial of state .with respect to 

Initial state elements 

2 H Partial of state with respect to 

Keplerlan state elements 

3 ■ Partial of state with respect to 

spherical state elements 

4 > Partial of state with respect to 

DODS state elements 



873-6660 

SPARES 




8 







NAMe, DESCRIPTION, ETC. 


TH: 

Tima (from epoch in seconds) of last 
acceleration in XDD 

H: 

Integrator stepslze, in sec&nds 

XDD<I,J): 

I-K-.10,K; J-1,3 

Satellite acceleration vectors 
K ■ acceleration array indicator 

TR£G(I): 

I-K-IO.K 

Time corresponding to acceleration 
array 

3X1(1); 

1-1,3 

First sum vector of satellite 
acceleration 

5X2(1): 

1-1,3 

Second sum vector of satelllLe 
acceleration 

XVDD(J1,J2,J3) ; J1-K-10,K; J2-1.3; J3-l,20 

Array of acceleraClon partlals 
K • accaleraclon array Indicator 

SV1(I,J): 

1-1,3; J-1,2U 

1st sum matrices for acceleration 
partlals 

SV2(I, J) ; 

1-1,3; J«l,20 

2nd sum matrices for acceleration 
partlals 

NSECTN 

Current section number 


ORIQIWAL PAGE tS 
OF POOR QUALITY 
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onmmi pass rs 

OF POOR QUALITY 


2.2 GPS Trajectory File (Unit 40) 


The GPS Trajectory File contains the positions and velocities (in 
ECI) of the GPS satellites. The number of satellites on the file is 
given on the header record and may vary from 1 to 24. 



RE 

BBli 



RECORD 

S 

wM 

BYTES 

NAME, DESCRIPTION, ETC. 


1 180 bytes 1-4 NB; Record Number - 1 


5-12 YEAR: Year of epoch of trajectory 

(2 digit) 

13-20 DAY: Day of epoch of trajectory 

21-28 SEC: Seconds of epoch of trajectory 

(UTC time) 

29-36 TVE-TO: Time of vernal equinox minus epoch 

of the trajectory 

37-44 ni: Interval at which trajectory is 

written (sec.) 

45-52 TLST; Last time on trajectory (sec. from 

epoch) 

53-60 STEP: Integration interval used when 

creating the trajectory (sec.) 

61-64 IFLOW: Kind of trajectory = 4, 5, 6 or 7 

(Not Used) 

65-68 150: Indicates if trajectory was made 

using epoch of date = 1 or epoch 
of 1950 = 0 

69-72 ICYC: Improvement cycle number (Not Used) 

73-80 DATE: Time clock value when trajectory 

was made (A8) 

81-84 NGPS: number of GPS satellites 

85-88 JGPS(l): Satellite number of 1st satellite 

(SV number) 

89-92 JGPS(2): Satellite number of 2nd satellite 


177-180 


JGPS(24): Satellite number of 24th satellite 
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OF took « 
POOR quality 







RECORD 

mu 

BYTES 

NAME, DESCRIPTION, ETC. 

All others 

212 

1-4 

NB: Consecutive record number = 2, 3, n 


to 

5-12 

TI : Seconds 

from epoch 


1500 





13-20 

X(l,l); X ) 

Inertial components of 



21-28 

X(2,l); Y 

satellite #1 position at 



29-36 

X(3,l); Z ) 

time TI (km) 



37-44 

XD(1,1): X 1 

Inertial components of 



45-52 

XD(2,1); Y 

satellite #1 velocity at 



53-60 

XD(3,1): Z ) 

time TI (km/sec.) 



61-68 

TB(1): Timing bias for satellite #1 (secsj 




(TBl 

= clock - truth time) 



69-76 

X(l,2); X ) 




77-84 

X 2,2 : Y } 

Satellite #2 position 



85-92 

X(3,2); Z ) 



93-100 

XD(1,2); X \ 




101-108 

XD(2,2): Y \ 

Satellite #2 velocity 



109-116 

XD(3,2): Z ) 




117-124 

TB(2): Timing bias for satellite #2 

• 




• 

XD(i,NGPS): 

^ ) 




XD(2,NGPS): 

Y > Satellite # NGPS velocity 




X0(3,NGPS): 

z ) 



12+56. NGPS 

TB(NGPS): Timing bias for satellite # NGPS 




ABCD(1,1): 

ABCD(2,1): 

1 Earth Centered Inertial 




• 

J. (19S0) to Earth Centered 




• 

i Fixed Transformation 




ABCD(3,3): 

) 




ABDCD(1,1): 

ABDCD(2,1); 

) Transformation of Velocities 




• 

> (XfTci * (ABDCD)^ . Xgp) 



156+56. NGPS 

ABDCD(3,3): 

) 
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Measurement Data File 


PASE FS 
OF POOR QUALITY 


The measurement data file input to PREFER is created by GTDS as the 
Observation Save File (GTDS Unit 46). This file contains all the informa- 
tion necessary to process the observations (e.g. station positions, refrac 
tion corrections, etc.). 


RECORD 


RECORD 

SIZE 



NSET: 


IPOINT (I,J): 


NAME, DESCRIPTION, ETC, 


Nuiaber of observation sets 
contained in data set 

I - 1,5; J - 1,7 

Pointer array of set related 

Information 


IPOINT (1,J) - Type indicator for set J 
J ■ 1, traclcing data 

• 2, telemetry data 

■ 3, PCE data 

• 4, relay data 

■ 3, landmark data 
■6 or 7, spare 



IPOINT (2,J) 
IPOINT (3,J) 
IPOINT (4,J) 

IPOINT (5,J) 


IPOINT a, a): 


miYPEd): 


■ pointer to reader record for 
set J 

* pointer to starting block for 
data from set J 

■ pointer to starting record 
within Indicated block for 
data from set J 

” Number of observations in set 
J 

I - 1,5 

Pointer array of set related 
Information (spare set) 

I - 1,30 

Array of external observation 
types to be used with input from 
this data set. Internal obser- 
vation types 1-30 correspond 
with IMTYPE(l) to IMTYPE(30). 
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or;giwal page is 

OF POOR QUALITY 


Ricono 



285-2S8 Sparc 
289-432 ISPRE U) : 


433-480 ISPRE (I) : 


481-576 RHDRl (I); 


577-720 RHDRl (I> ; 


721-864 RHDR2 (I) ; 


865-960 RHDR2 (I) : 


961-1008 RHDR3 (I) : 


1009-1152 RHDR3 (I): 


1153-1200 RHDR3 (I) 


1201-1296 RHDR4 (I) : 


1297-1440 RiroR4 (I): 


NAMI, DSSCRimON. BTC 


I - 2,37 
Spare locations 

I - 38,49 
Spare locations 

I - 1,24 

Header Information for first 
observation sec 

I - 25,420 

Header Information for first 
observation sec 

I - 1,396 

Header Infors-itlon for second 
observation set 

I.- 397,420 

Header Information for second 
observation set 

I - 1,12 

Header Information for third 
observation sec 

I - 13,408 

Header Information for third 
observation sec 

I - 409,420 

Header Information for third 
observation sec 

I - 1,24 

Header information for fourth 
observation set 

I - 25,420 

Header Information for fourth 
observation set 


1441-7200 Spare 
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GF p 


NAMf . OISCRIPTION, ITC. 


RHDR5 

(I): 

I - 1,420 


RHDR6 

(I): 

I - 1,420 


RBSR7 

(1)5 

I - 1,420 


RHDR8 

(I); 

I - 1,396 

Headers for observation sets 5 
through b 

!' 

! j 

I ' 

RHDR8 

(I): 

I - 397,420 

Header information for observe'* 
tlon set 8 (continued) 

t’ 

I 

? 

il 


Spar* 


For a dascrlptlon of header Information for an 
obaervaClon set« see data set layout for GTDS 
Observation Tape* FORTRAN logical unit 29[4] . 



*NOTE ; Record 2 exists in this format 

only If there are more than four 
observation sets on the file. 

If there are four or less, record 
2 Is formatted like records 3-N 
(data records) . 


f 

i' 


' fS' 

7Y 
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Of POOH QUALIVY 


RECORD 

RECORD 

SIZE 

BYTES 

NAME, DESCRIPTION, ETC. 

ALL 

(data) 

7200 

1-8 

TOBYMD; 

Corrected observation time tag 
(YYMMDD. ) 



9-16 

TOBHMS: 

Corrected observation time tag 
(HHMMSS.SSS) 



17-24 

OMl: 

Uncorrected observation 
measurement No. 1 



25-32 


Doppler Count Interval (sec.) 



33-40 

SIGl: 

Standard deviation of data 
reduction fit 



41-48 

OMl A: 

Corrected observation measurement 
No. 1 



49-56 

IBIT{8): 

Validity flags 



57-64 

DMR: 

Observation refraction correction 



65-72 

TOBSC: 

TOBYMD and TOBHMS converted to A.l 
seconds from Jan. 1,; DC closes 
file zero hour of epoch year 



73-80 

DAMTD: 

Sum of antenna mount and 
transponder delay correction 



81-88 

JBIT(8): 

Identi fiction flags 



89-96 

TCOR: 

or 

TOBSC corrected for bias and light 
time (in ET seconds from epoch). 




FREF: 

For TORS data, primary frequency 
(MHz). 



97-100 

NOBS: 

Working file observation number 



101-104 

IGATE; 

Range-Gating indicator 



105-108 

MTYPE: 

Observation type Indicator 


15 
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ORIGINAL PAGE IS 

“ I Al 


RECORD 

RECORD 

SIZE 

BYTES 

NAME, DESCRIPTION, ETC. 

ALL 

(data) 

7200 

105-108* 
(TDRS Data 
only) 

MTYPE; 

Packed indicator 
■ uplink ID * 1000000. + ground 
transponder number * 10000. 

+ downlink ID * 100. + MTYPE 
(ID is 2-digit referenced to 
60 byte format) 



109-110 

I SOU: 

Observation source indicator 



111 

IBF: 

Doppler Bias Frequency indicator 
for SST data. Not used for TDRS 
data. 



112 

IPFO: 

Primary Frequency Offset indicator 
for SST data. Not used for TDRS 
data. 



113-116 

lEDIT: 

Edit flag 



117-120 

I BIAS; 

Bias type indicator 



121-124 

ISTA; 

Station geodetic indicator 



125-128 

lUPDTE: 

or 

Record update indicator 




100*ISTAAT 
+ISTAA: (TDRS only) 

ISTAAT * Ground transponder 
internal index 
(1,2,3...) 

ISTAA = Transmit station internal 
index (1, 2, 3 ...) 



129-132 

IFGP(I); 

1*1,4; FLags indicating 
corrections made prior to GTDS/ODS 
processing 



133-136 

WORD10: 

WORD 10 of 60 byte record 



137-140 

IFGB(I); 

1*1,4; Flags indicating 
corrections made during GTDS/ODS 

processing 



141-144 

W0RD9: 

Word 9 of 60 byte record 



145-7200 

49 additional observations containing 
duplication of bytes 1-144 definition 
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2.4 Solar/Lunar/Planetary (SLP) (Unit 20 and 22) 

There are two SLP files available to PREFER on the GSFC 360/95. One 
uses mean of 1950.0 coordinates (Unit 20) while the other uses true-of-date 
coordinates (Unit 22). PREFER automatically selects the correct file to 
match the satellite ORBIT file. The current namas of the two SLP files 
are: 0RBIT.GTDS.SLP1950.DATA and ORBIT.GTDS.SLPTOO.DATA. 

2.5 Time Coefficients File (Unit 21) 

The time coefficients file contains Information on the time system 
transformations and polar motion. The current name of the file on the 
GSFC 360/95 Is ORBIT.GTDS.TIMCOF.DATA. 

2.6 TORS ORBIT Files (Units 41. 42, 43) 

When TORS relay measurements are processed, the TORS ephemerldes are 
obtained from three, direct access, ORBIT files on the GSFC 360/95. The 
format of these files Is Identical to that of the Satellite ORBIT File 
(Unit 50) except for the difference between sequential and direct access. 
There 1s no unique relationship between the unit number (41, 42 or 43) and 
the Individual TORS: PREFER determines this relationship from the 7-d1g1t 
satellite Identifier on the header of the TORS ORBIT files and on the TDRID 
cards. 


It Is Important that the TORS ephemeris coordinate system (1950 or 
true-of-date) match the coordinate syst'jm on the Satellite ORBIT file. If 
they do not match, the PREFER job will terminate. 

PREFER does not actually use the TORS ORBIT files when processing 
measurements. At the beginning of the job, PREFER Interpolates the ORBIT 
files (at a rate determined by the STEP card) and copies the ephemerldes to 
a file (also on unit 40) with the same format as a GPS Trajectory file. 

This Is the file used to process the measurements. 
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3.0 CARD INPUT 


All card input to PREFER consists of 72-column keyword cards which 
have seven variables per card and are read using IBM's free form input 
(except for the keyword): 

Variable Type 

1 Alphanumeric keyword (this must be located in columns 
1-8, left justified) 

2 Integer 

3 Integer 

4 Real 

5 Real 

6 Real 

7 Real 

The last six variables may be located anywhere in columns 9-72 
since blanks or commas are used to separate the fields. The order of 
the variables is important and thus unused variables must still be input 
(a zero or a double comma is suggested). See the IBM Fortran (extended) 
manual for further information. 

There are two general classifications for card types. Run control 
cards are used to input constants required for the run. Run initialization 
cards are used to set a pviovi standard deviations, state noise standard 
deviations and time constants of Markov processes. These cards also deter- 
mine the parameters to be included in the state vector. 

With the exception of the STEP card, any field which is blank will be 
read as a zero and will override the default values (if any). There are 
no restrictions on the order in which cards are read. 

Before describing the format of the individual cards, it is necessary 
to discuss the dynamic models and the requirements for the card input. 

PRECEDING. PAGE BLANK NOT, FILMED’ 
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PREFER has the capability to solve for the following set of parameters 

satellite x-position at epoch 

2 satellite y-position at epoch 

3 satellite z-position at epoch 

4 satellite x-velocity at epoch 

5 satellite y-velocity at epoch 

6 satellite z-velocity at epoch 

7 satellite drag coefficient 

8 perturbing gravitational acceleration (vertical) 

9 perturbing gravitational acceleration (cross-track) 

10 perturbing gravitational acceleration (along- track) 

11 acceleration of 1st thrust segment (vertical) 

12 acceleration of 1st thrust segment (cross-track) 

13 acceleration of 1st thrust segment (along-track) 

14 acceleration of 2nd thrust segment (vertical) 

15 acceleration of 2nd thrust segment (cross-track) 

16 acceleration of 2nd thrust segment (along-track) 

17 host satellite clock timing error 

18 host satellite clock drift rate 

19 altimeter bias 

20 ground station measurement bias 1 

21 ground station measurement bias 2 

22 ground station refraction parameter 

23 ground station position error (x) 

24 ground station position error (y) 

25 ground station position error (z) 

26 GPS satellite position error (H) 

27 GPS satellite position error (C) 

28 GPS satellite position error (L) 

29 GPS satellite clock timing error 
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All the parameters (except the first six) are optional. That is, 
the user may elect to solve for almost any subset of these parameters in 
a given run. The only restrictions are that parameters which appear in 
groups of three (e.g., gravitation, thrust, ground station positions and 
GPS positions) will be treated as a group (all three will be adjusted or 
none will be adjusted). 

The first nineteen of these parameters are called dynamic parameters 
because their values may change with time. This may happen simply because 
they have a non-identity transition matrix or because the parameters are 
driven by some random process (i.e., they have state noise). Of these 
nineteen parameters, six (drag, gravitation, clock drift and altimeter 
bias) are assumed to be first order Markov processes; i.e., if X represents 
one of these parameters, then 


1 *. 

T 


where 


T = the correlation time 
u = zero mean white noise . 

For these six parameters, it is necessary for the user to input three 
separate quantities: the standard deviation of the error in the a pHoH 

parameter estimate, the standard deviation of the process output (a^) and 
the correlation time. For all the other parameters, (inly the a pviovi 
standard deviation and the state noise spectral density are required. 

Determination of the appropriate values for the a priori sigma and 
the correlation time is usually not difficult. However, calculation of 
the state noise is frequently done with a bit of hand waving. In fact, 
the bias that some people have against Kalman filtering is often the result 
of their misunderstanding of the meaning of state noise (and, consequently, 
their inappropriate choice of values). This is somewhat understandable 
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since few books elaborate on this subject. The primary purpose of state 
noise is to account for errors in the dynamic model. 

In the PREFER program, the spectral density of u is computed internally 
from the state noise sigma that the user inputs. That is, the user specifies 
the standard deviation of x (in the previous equation). The standard 
deviation on x that results from integration of the first-order Markov 
process is related to by the equation 



Therefore, will be computed internally by the program. 

For the orbital elements and thrust accelerations, (the non-Markov 
processes), the input quantity for state noise is the expected standard 
deviation of the time derivative of the parameter. In this case, the standard 
deviation of the parameter will grow as where t is the integration 

time. Since the thrusting times will generally be short, it should not 
be necessary to input any state noise for thrust. However, the use of 
a very small value may help minimize numerical problems. 

A similar situation also exists for the orbital elements. If no other 
force model parameters were adjusted, then it would be necessary to include 
state noise on the orbital elements (particularly the velocity terms). 

However, if the gravitation and/or other parameters are adjusted, then 
the state noise on the orbital elements should just be large enough to 
alleviate numerical problems. 

Notice that no provision was made for specifying state noise on the 
host clock timing error. This was done because the error is really in 
the clock frequency (drift rate) and, thus, errors in the clock time can 
be obtained exactly by integrating the drift rate. 

All the remaining parameters are related to the measurements and are 
not part of the dynamic model. They are all assumed to be constant (identity 
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state transition matrix) and also have no state noise. For these reasons, 
they come under the general heading of ^'measurement biases". 

These parameters are treated differently than dynamic parameters In 
several aspects but the most significant difference Is their use as "pass 
parameters". That Is, parameters for a particular station or satellite do 
not actually enter the filter until measurements from that station 
(satellite) are processed. Once the pass Is over, the parameters are 
dropped freom the filter. This procedure has been shown to produce the 
same results for the common parameters as If all parameters were carried 
throughout the run (assuming that all pass parameters are Independent 
between passes). 

In the case of station position errors, pass Independence Is not a 
valid assumption. Thus, the assumption of pass Independence means that 
some Information Is not being used and the filter Is slightly suboptimal. 
However, the filter estimate will not be biased because of this assumption. 

It should be noted that the original design of PREFER did not attempt 
to solve for station and GPS position errors and GPS timing errors; they 
were carried as "consider" parameters. However, problems In Implementing a 
smoother with consider parameters forced their Inclusion as adjusted 
parameters. At this time, It 1s not known how severe the numerical 
problems are going to be when many pass parameters are simultaneously 
adjusted. It was assumed that a maximum four ground satlons and 15 GPS 
satellites* could be simultaneously observable. If all possible parameters 
were adjusted, this would result In a state dimension of ^ (PREFFR allows 
for 100 states). Numerical problems should not be unexpected with a state 
vector of this size. 

The user should also be wary of runs which simultaneously adjust the 
drag coefficient and gravitational accelerations. If the orbit Is nearly 
circular, the drag acceleration will probably not be distinguishable from 
the along-track gravitational acceleration (unless the correlation times 
are significantly different). 

* PREFER also has a total limit of 50 ground stations and 24 GPS satellites 
In a run. 
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The description of the STEP card (Section 3.2.10) uses the expression 
"mini -batch measurement processing". In order to reduce the program 
running time and disk storage, measurements are processed in small batches, 
during which time (typically 120 seconds), the dynamic model errors are 
assumed to be negligible. State noise is included only when the filter 
proceeds from one mini -batch to the next. 

The processing of TORS relay measurements is handled somewhat 
differently than that for GPSPAC measurements. First, it is not possible 
to estimate the TORS ephemeris errors using the current version of PREFER. 
Secondly, it is necessary for the user to establish a unique relationship 
(using the TDRID card) between the TORS and ground tracking stations. It 
is assumed that each of these ground stations tracks only one TDRS and 
nothing else. Thus the deletion of the station position errors and 
measurement biases (pass parameters) from the state vector is determined by 
the visibility of the user satellite as viewed from the TORS (see the ELCllT 
card description for details). 
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3.2 Run Control Cards 


All the following cards are optional. 


Keyword Description 


3.2.1 

DEBUG 

Turns on debug print 

3.2.2 

DRA6C0N 

Coefficients used for atmospheric density modelling 

3.2.3 

EARTH 

Earth ellipsoid parameters, and speed of light 

3.2.4 

ELCUT 

Measurement elevation cutoffs 

3.2.5 

MED IT 

Editing threshholds 

3.2.6 

MSIG 

Overrides measurement sigmas on measurement file 

3.2.7 

PRINT 

Printout options 

3.2.8 

SMR 

Satellite area to mass ratio 

3.2.9 

STAPOS 

Ground station positions 

3.2.10 

STEP 

Mini-batch step size, stop time and integration step 
size 

3.2.11 

TDRID 

Defines the relationships between TDRS Identifiers and 
ground tracking stations 

3.2.12 

THRTIM 

Times of thrust events 


3.2.1 DEBUG 


original: page is 

OF POOR QUALITY 


The DEBUG card turns on the debug printout for specified subroutines. 
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Variable Digit Type Description 

1 Alphanumeric (A8) Keyword DEBUG. 

2 1 7 digit integer 1 in this digit will turn on the debug 

print for subroutine FILTER. 

2 1 in this digit will turn on the debug 

print for subroutine MEAS. 

3 1 In this digit will turn on the debug 

print for subroutine CARD. 

4 1 In this digit will turn on the debug 

print for subroutine DYNAM. 

5 1 in this digit will turn on the debug 

print for subroutine SNOISE. 

6 1 In this digit will turn on the debug 

print for subroutines GETHDR and INTERG. 

7 1 In this digit will turn on the debug 

print for subroutine EVAL and EARTH. 
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DEBUG 


Notes: (1) The default Is no debug print. 

(2) A brief description of the debug printout is given in the 
Appendix. In general, it would be necessary to have a 
Fortran listing of PREFER to interpret the output. 
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ORIGINAL' PAGE B 

3 • 2 • 2 DRAGCON OF POOR QUALITY 

The DRAGCON card allows the user to input the five coefficients which 
are used to model the density of the earth's atmosphere. 

IM6GQrrtl08«3 0 ,01962 -«,93S9 , 0001018 
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Variable Digit Type 


Description 


1 

2 

3 

4 

5 

6 


1-2 

3-4 

5-6 


Alphanumeric (A8) Keyword DRAGCON. 

Integer Index (i) of first coefficient on card. 

Integer Index (j) of second coefficient on card. 

Integer Index (k) of third coefficient on card. 

Integer Not used. 

Real Coefficient corresponding to number in 

digits 1-2 (d.) . 

Real Coefficient corresponding to number in 

digits 3-4 (d.) . 

J 

Real Coefficient corresponding to number in 

digits 5-6 (dj^) . 


Note: The five coefficients are used in the equation 


p ■ exp(d.|h - d 2 - i^^h^ + d^h - dg ) 

3 

where h is the altitude in kilometers and p is density in kg/km . 

The defaults for d-| • dg are shown on the sample card. The numbers in 
columns 9>10 and 11*12 should not be zero (or blank). 
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3.2.3 EARTH 


The EARTH card is used to input the Earth's semi -major axis, 
eccentricity, gravitational constant and the speed of light. 
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Variable Type 


Description 


1 

2 

3 

4 

5 


6 

7 


Alphanumeric (A8) Keyword EARTH. 

Integer Polar motion switch: 0 = no polar 

motion applied, 1 = polar motion 
applied. Default = 1. 

Integer Not used. 

Real Earth semi -major axis (km). 

Default = 6378.166. 

Real Earth eccentricity or inverse of 

flattening. If value is greater than 
1, program assumes 1/f. 

Default = .081813334(e) or 298,3(l/f) 

Real Earth gravitational constant (GM) in 

km^/sec^. Default = 3.9860064x10^. 

Real Speed of light (km/sec). 

Default = 2.997924562x10®. 
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3.2.4 ELCUT ORIGim PAGE rs 

OF POOR QUALITY 

The ELCUT card 1s used to input the elevation cutoffs for processing 
ground measurements and GPS measurements. 
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Variable 

Type 

Description 

1 

Alphanumeric (A8) 

Keyword ELCUT 

2 

Integer 

Not used. 

3 

Integer 

Not used. 

4 

Real 

Elevation cutoff (degrees) for processing 
ground tracking measurements^. All 
measurements below this cutoff will be 
deleted. Default = 5“. 

5 

Real 

Elevation cutoff (degrees) as measured 
at the host satellite for processing 
GPS measurements. All measurements 
below this cutoff will be edited. 
Default =0®. 
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ELCUT 


Variable Type Description 

6 Real Radius of sphere surrounding earth which is used to 

determine whether a user satellite is "visible” to a 
TORS (See note 2). 

Default = 6778.0 km. This visibility test determines 
when pass parameters associated with TDRS tracking 
(e.g. station biases, position errors, etc.) are 
deleted from the filter state. 

Notes: 1) The formula for computing refraction effects will blow up if 

elevation is less than or equal to zero. Therefore, the 
elevation cut-off should be greater than zero. 

2) The shaded area in the figure is not "visible" from the TORS. 


Visibility radius 
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mil 


3.2.5 MEDIT 


The MEDIT card is used to set the editing and printing thresholds 
for outlying measurements. 
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Variable 

Type 

Description 

1 

A6 

Keyword MEDIT 

2 

Integer 

Not used 

3 

Integer 

Not used 

4 

Real 

Measurement editing threshold. Default = 20 
Measurement is also printed. 

5 

Real 

Measurement printing threshold. Default = 10 
Measurement is printed but not edited. 


NOTE: (1) Editing or printing of outlying measurements will occur when 

2 2 9 

the ratio z /a^ (where z is the measurement residual and 

is the residual variance computed by the filter) exceeds the threshold. 

Note that unlike most batch processors, the procedure includes the 

uncertainty in the state estimate as well as the measurement noise. 
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MSIG 


3.2.6 MSIG 


The MSIG card 1s used to override the irieasurement standard deviations 
(as given on the measurement tape) for particular stations* GPS satellites 
or satellite altimetry. 
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Variable Type 


Descri pti on 


1 

2 

3 


4 


Alphanumeric (A8) Keyword MSIG. 

Integer Station or satellite number (maximum of 

5 digits and must be positive). 

Integer Source of measurement. 

1 = GPS satellite 

2 = Ground station 

3 = Altimeter 

Real Measurement standard deviation (kilometers 

or dimensionless) for the fir^ measure- 
ment from a station, satellite or alti- 
meter. Default = SIGl on Measurement 
Data File (see page 15). 
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Variable Type Description 

5 Real Measurement standard deviation (kilometer or km/sec) for 
the second measurement from a satellite which has a pair 
of measurements. This only applies to the psuedo-delta 
range from GPS satellites and range rate for ground 
stations (including TORS relay measurements). 

Default ■ SI61 on Measurement Data File. 


Note: Measurements processed in pairs are: 


Ground range and range rate: 

GPSPAC 

TORS 


Measurement 1; 
range 

psuedo range 
range 


Measurement 2 
range rate 
psuedo delta- range 
range rate 


i;¥.©E B 

OF FOQR QUAlI fY 
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3.2.7 PRINT ORIGIMAL PAGE IS 

OF POOR QUALITY 

The PRINT card defines the printout options. 
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Variable 

1 

2 


Digit Type Description 

Alphanumeric (A8) Keyword PRINT. 

Integer 5 digit integer. 

1 Defines coordinate system of printed and 

plotted output: 0 = HCL, 1 » i-CI. 

Default = 0. 

2 Controls filter print: D = no state 

output, 1 = orbital elements only, 

2 » full printout. Default = 2. 

3 Controls smoother print: 0 = no state 

output, 1 = orbital elements only, 

2 = full printout. Default = 2. 

4 Controls generation of the smoother 

covariance (EC I coordinates) on unit 71. 

0 » no output, 1 = output. Default = 0. 
NOTE: The computer I/O charges will 

approximately double when using 
this option. 
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Variable Digit Type Description 

r 

2 5 Integer Controls plotting of filter/ smoother 

standard deviations on position. 

0 ■ no plots, 1 » filter plot only, 

2 - smoother plot only, 3 * filter and . 

smoother plots. Default ■ 0. 
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OF POOR QUALITY 


SMR 


3.2.8 SMR 

The SMR card 1s used to input the satellite area to mass ratio. 
4MR 0 0 0,311E-8 ,,,, 

t <r 
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Variable Type Description 

1 Alphanumeric (A8) Keyword SMR. 

2 Integer Not used. 

3 Integer Not used. 

4 Real The satellite area to mass ratio in 

2 

kilometers /kilogram. Default is the 
ratio of area to mass on the ORBIT 
file. 


38 


Bvsii^sand TECjmpuiGiC4L Systems, Inc 


STAPOS 


'?.1G5WAU PAQE 5S 

3.2.9 STAPOS POOR QUALITY 

The STAPOS card is used to input the positions of the ground tracking 
stations. 
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Variable Type 


Description 


1 

2 

3 

4 

5 

6 


Alphanumeric (A8) 

Integer 

Integer 

Real 

Real 

Real 


Keyword STAPOS. 

Station number (maximum of 5 digits). 

Not used. 

Station geodetic latitude (degrees). (|) 

Station East longitude (degrees). X 
Station height above ellipsoid (kilometers), h 


NOTE: (1) It should not be necessary to input any station positions using 

the STAPOS card; the header record on the Observation Save File 
(Unit 9) contains all the station positions. However, if the 
user wanted to intentionally change station positions, this can 
' be done by using the STAPOS card. 
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STEP 


3.2.10 STEP ORm'Al PAGE IS 

OF POOR QUALITY 

The STEP card controls the step size used for filter/ smoother output 
and also controls the stop time of the run. 

•STEP 0 I) 120,0 86400,0 120, 600, 

# • 
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Variable Type 


Description 


1 A6 

2 Integer 

3 Integer 

4 Real 


5 Real 


Keyword STEP. 

Not used. 

Not used. 

Step size (seconds) for mini -batch 
measurement processing. All filter/ 
smoother output will use this interval (see note 1) 
Default is 120 sec. 

Stop time (seconds) from epoch of host 
trajectory tape. The default is 86400 
seconds but the program will stop at the 
end of either the host trajectory tape or 
the measurement tape. 
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STEP 


Variable Type nescription 

6 Real Maxitnum step size for use In the Taylor series 

Integrator. 

Default = 120 sec. 

7 Real Step size used In copying TDRS ORBIT files to a 

TDRS trajectory file. 

Default » 600 sec. 

Notes: (1) This step size should be chosen so that dynamic modelling 

errors are negligible during the interval because the filter 
assumes that state noise is zero. 

(2) If zeroes or blanks are input on this card, the program will 
use the default value. 

(3) There is an additional restriction that the stop time divided 
by the mini-batch step size Is less than 999. This restric- 
tion is imposed because PREFER must have sufficient core storage 
to store the pointers used in the random-access disk I/O. 
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ORIGINAL’ PAGE IS 

3.2.11 TDRID OF POOR QUALITY 

TDRID cards are used to define the relationships between the 
seven-digit and two-digit TORS identifiers and ground tracking stations. 

These cards must be input when TORS relay tracking is processed. 

TDRID 37 II 3, XXXJWX, »»» 
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Variable Type 

1 A5 

2 Integer 


3 Integer 

4 Real 


5 Real 


Description 

Keyword TDRID. 

Tracking station number (two-digit). 
This number must match the station 
number on the header of the Measurement 
Data File. 

Not used. 

Two-digit TDRS identifier. This number 
must match the TDRS number on the 
Measurement Data File. 

Seven-digit TDRS identifier. This 
number must match the TORS number on 
the TDRS ORBIT file. 
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3.2.12 THRTIM 


ORIGINAL’ PAGE FS 
OF POOR QUALITY 


The THRTIM card is used to input the start and stop times for two 
thrusting periods. If omitted, no thrusting will be used. 
fHRTIW 0 0 5700, 5760, 10500, 10590, 
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Variable 


Jm 


Description 


1 A6 

2 Integer 

3 Integer 

4 Real 

5 Real 

6 Real 

7 Real 


Keyword THRTIM. 

Not used. 

Not used. 

Start time (seconds from epoch) 
of first thrust segment. 

Stop time (seconds from epoch) of 
first thrust segment. 

Start time (seconds from epoch) of 
second thrust segment. Default = 0. 
Stop time (seconds from epoch) of 
second thrust segment. Default = 0. 
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3.3 Run Initialization Cards Op POOR QUALITY 

All the following cards are optional but it may be unreasonable to 
make runs without initializing certain parameters. 

Because the cards may be read in any order, it is important that the 
CLOCK and CLOCKQ cards have the same number for variable 3. 



Keyword 

Description 

3.3.1 

ALTIM 

a priori a for altimeter bias 

3.3.2 

ALTIMQ 

state noise for altimeter bias 

3.3.3 

ALTIMT 

de-correlation time for altimeter bias 

3.3.4 

CLOCK 

a priori a for host clock timing bias and 
clock drift rate 

3.3.5 

CLOCKQ 

state noise for host clock timing bias and clock 
drift rate 

3.3.6 

CLOCK! 

de-correlation time for clock drift rate 

3.3.7 

CSTA 

a priori a for station position errors 

3.3.8 

DRAG 

a priori a for drag coefficient 

3.3.9 

DRAGQ 

state noise for drag coefficient 

3.3.10 

DRAG! 

de-correlation time for drag coefficient 

3.3.11 

GPSP 

a priori a for GPS position errors 

3.3.12 

GPS! 

a priori a for GPS timing errors 

3.3.13 

GRAV 

a priori a for perturbing gravitational 
accelerations 

3.3.14 

GRAVQ 

state noise for perturbing gravitational 
accelerations 

3.3.15 

GRAVT 

de-correlation time for perturbing gravitational 
accelerations 

3.3.16 

MBIAS 

a priori a for station biases 

3.3.17 

REFRAC 

a priori a for refraction parameters 

3.3.18 

STATEl 

a priori a for first three orbital elements 

3-3.19 

STATE! Q 

state noise for first three orbital elements 

3.3.20 

STATE2 

a priori a for orbital elements 4-6 

3.3.21 

STATE2Q 

state noise for orbital elements 4-6 
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Keyword 

Description 

3.3.22 

THRUSTl 

a pvioTyi a 
segment 

3.3.23 

THRUST! Q 

state noise 
segment 

3.3.24 

THRUST2 

a pviorn, a 
segment 

3.3.25 

THRUST2Q 

state noise 
segment. 


for acceleration of first thrust 
for acceleration of first thrust 
for acceleration of second thrust 
for acceleration of second thrust 


Because certain parameters require more than one input card (e.g., 

GRAV, GRAVQ, GRAVT), it is possible to have unusual errors occur if the 
deck is not correctly set up. For example, suppose that the GRAVQ card was 
input but no GRAV card was provided. The presence of the GRAVQ card tells 
PREFER that gravitational parameters are to be adjusted. However, the 
a prioH variances will default to zero because no GRAV card was supplied. 
This will obviously result in problems for the filter. 
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3.3.1 MM ORiGIMAL PAGE IS 

OF POOR QUALITY 

The ALTIM card is used to initialize the standard deviation for the 
altimeter bias. If omitted, the altimeter bias is not adjusted. 

ALTin 00 0,005 ,,,, 
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Variable Type 


Descri ption 


1 Alphanumeric (A8) Keyword ALTIM. 

2 Integer Not used. 

3 Integer Not used. 

4 Real A priori standard deviation for the 

altimeter bias (km). No default. 
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ALTIMQ 


ORIGINAL PAGE \B 

3.3.2 ALTIMQ OF POOR QUALITY 

The ALTIMQ card is used to Input the state noise on the altimeter 

bias. 
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Variable Type 


Description 


1 

2 

3 

4 


Alphanumeric (A8) Keyword ALTIMQ. 

Integer Not used. 

Integer Not used. 

Real Steady sigma of the altimeter bias 

Markov process (in km). Default = 0. 


48 


BVSINFSiANO TECMmumCAL SYSTEMS, INC 


3.3.3 ALTIMT 


omm P/ice ,s 

O' POOR QU/iuTy 


The ALTIMT card is used to input the time constant of the first-order 
Markov process for the altimeter bias. 
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Variable 


Descri ption 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Keyword ALTIMT. 

Not used. 

Not used. 

Correlation time (seconds) for changes 
in the altimeter bias. Default = 1.0 
seconds . 
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CLOCK 


0RK3INAU PAGE fS 

3.3.4 CLOCK QUALiTY 


The CLOCK card is used to Initialize the standard deviations for the 
host clock timing error and drift rate. It omitted, the v'.lock errors 
will not be adjusted. 
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Variable Type 


Description 


1 

2 

3 


4 


5 


Alphanumeric (A8) 

Integer 

Integer 


Real 


Real 


Keyword CLOCK. 

Not used. 

Number of parameters: 

1 * adjust clock drift rate only 

(e.g. NAVPAC). 

2 = adjust clock timing error plus drift rate 

(e.g. GPSPAC rate). 

A priori standard deviation for host clock 
drift rate (milliseconds/second). No 
default. 

A priori standard deviation for host clock 
timing bias (milliseconds). No default, but 
this number is not used if column 15 is 1. 
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origim/^l page (S 

OF POOR QUALITY 


CLOCKQ 


3.3.5 


CLOCKQ 


The -KQ card Is used to Input the state noise for the clock drift 

rate. 
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Variable Type 


Description 


1 Alphanumeric (A8) 

2 Integer 

3 Integer 


4 Real 


Keyword CLOCKQ. 

Not used. 

2 

Number of parameters j 

1 * adjust clock drift rate only 

2 = adjust clock timing error plus 

drift rate. 

Steady state sigma of the host clock 
drift rate, Markov process^ (in mi Hi 
seconds/sec). Default * 0. 


Note: (1) The Markov process only applies to the clock drift rate . The 
clock timing error is obtained by integrating the rate. 

(2) This number should match the number on the CLOCK card. 
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CLQCKT 


3.3.6 OOCKl 

The CLOCKT card Is used to Input the time constant of the first order 
Markov process for the host clock drift rate. 

OLDCKT 0 2 50000, »>»* 
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Variable Type 


Description 


1 

2 

3 


4 


Alphanumeric (AS) 

Integer 

Integer 


Real 


Keyword CLOCKT. 

Not used. 

Number of parameters: 

T = adjust clock drift rate only 
2 = adjust clock timing error plus 
drift rate. 

Correlation time (seconds) for changes 
in the host clock drift rate. Default - 
1 .0 seconds. 


Note: The Markov process only applies to the clock drift rate . The 
clock timing error is obtained by integrating the rate. 
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CSTA 


3.3.7 CSTA ■ OmClWAL PAGE IS 

OF POOR QUALITY 

The CSTA card Is used to Initialize the standard deviation for the 
ground station position errors (In earth-centered fixed Carteslon 
coordinates). If omitted, ground station position errors will not be 
adjusted. 
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Variable Type 


Description 


1 Alphanumeric (A8) 

2 Integer 

3 Integer 

4 Real 


5 Real 


6 Real 


Keyword CSTA. 

Station number (if 0, this card will 
apply to all stations). 

Not used. 

A priori standard deviation for the X 
component of station position error (km). 
No default. 

A priori standard deviation for the y 
component of station position error (km). 
No default. 

A priori standard deviation for the z 
component of station position error (km). 
No default. 
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CSTA 


Note: (1) If it is desired to use the same sigmas for all stations but 

one or two* then use one card with no station number and the 
other cards with station numbers. 
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ORlOilMAL PAOi II 
OF POOR QUALITY 


DRAG 


3.3.8 DRAG 


The DRAG card is used to initialize the standard deviations for the 
drag coefficient. If omitted, the drag coefficient is not adjusted. 
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Variable Type Description 

1 Alphanumeric (A8) Keyword DRAG. 

2 Integer Not used. 

3 Integer Not used. 

4 Real A priori standard deviation for the 

satellite drag coefficient (dimensionless). 
No default. 
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OF POOR QUALITY 


DRAGQ 


3.3.9 DRAGQ 


The DRAGQ card is used to input the state noise for the drag coefficient. 
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Vari abl e Type Description 

1 Alphanumeric (A8) Keyword DRAGQ. 

2 Integer Not used. 

3 Integer Not used. 

4 Real Steady state sigma of the drag 

coefficient Markov process (dimension- 
less). Default = 0. 
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DRAG! 


3.3.10 DRAG! POOR QUALITY 

The DRAGT card is used to input the time constant of the first-order 
Markov process for the drag coefficient. 
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Variable 


1 
2 

3 

4 

in the drag coefficient. Default = 
1 .0 seconds. 


Type Description 

Alphanumeric (A8) Keyword DRAGT. 

Integer Not used. 

Integer Not used. 

Real Correlation time (seconds) for changes 
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GPSP 


3.3.11 GPSP 


The GPSP card is used to initialize the standard deviations for the 
GPS position errors (in earth-centered local corrdi nates). If omitted, 
GPS position errors will not be adjusted. 
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Vari abl e Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 
Integer 

Integer 

Real 


Real 


Real 


Keyword GPSP. 

GPS satellite number (if 0, this card will 
apply to all satellites). 

Not used. 

A pHoz^ standard deviation for the 
H (vertical) component of GPS position 
error (km). No default. 

A priori standard deviation for the C 
(crosstrack) component of GPS position 
error (km). No default. 

A priori standard deviation for the L 
(alongtrack) component of GPS position 
error (km). No default. 


Note: (1) Note (1) on the CSTA card also app 

(2) GPS position errors are treated as 
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ORIGIWAL PAGE IS 
OF POOR QUALITY 


GPST 



3.3.12 GPST OKIGIIMAL PAGE IS 

OF POOR QUALITY 

The GPST card is used to initialize the standard deviations for GPS 
timing error. If omitted, GPS timing errors will not be adjusted. 


GPST 


a 0 


0, 000003 
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I Variable Type 

1 Alphanumeric (A8) 

2 Integer 

W 

3 Integer 

I ! 4 Real 

Note: 


Description 

Keyv/ord GPST. 

GPS satellite number (if 0, this card 
will apply to all satellites). 

Not used. 

A priori standard deviation for the GPS 
timing error (milliseconds). No default. 


(1) Note (1) on the CSTA card also applies. 

(2) GPS position errors are treated as "pass" parameters. 
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6RAV 


3.3.13 6RAV 


The 6RAV card is used to initialize the standard deviations for the 
perturbing gravitational accelerations. If omitted, the perturbing 
gravitational accelerations are not adjusted. 


oPAV 0 0 S. 4 E -9 0 , 7 E -9 l.OE -9 m 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword GRAV. 

Not used. 

Not used. 

A priori standard deviation for the 
vertical component of gravitational 
acceleration (km/sec ). No default. 

A priori standard deviation for the cross- 
track (perpendicular to the orbit plane) 
component of gravitational acceleration 
(km/sec ). No default. 

A priori Standard deviation for the along 
track component of gravitational accelera- 
tion (km/sec ). No default. 
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GRAVQ 


3.3.14 GRAVQ 


The GRAVQ card is used to input the state noise for the perturbing 
gravity accelerations. 


6PW0 0 0 S.4E-9 0.7E-9 l.OE-9 m 
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Van" abl e 


Type 


Description 


1 

2 

3 

4 


Alphanumeric (AS) 

Integer 

Integer 

Real 


Real 


Real 


Keyword GRAVQ. 

Not used. 

Not used. 

Steady state sigmas of the vertical 
gravitational acceleration Markov 
process (km/sec^). Default = 0. 
Steady state sigmas of the cross - 
track gravitational acceleration 
Markov process (km/sec**). Default = 
Steady state of the along-track 
gravitational acceleration Markov 
process (km/sec^). Default = 0. 
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3.3.15 GRAVT 


OF POOR QUAUTY 


The GRAVT card is used to input the time constants for the first-order 
Markov processes of the perturbing gravitational accelerations. 


6PAVT 0 0 iPOO, POO, POO, M 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) Keyv/ord GRAVT. 

Integer Not used. 

Integer Not used. 

Real Correlation time (seconds) for changes 

in the vertical gravitational accelera- 
tions. Default = 1.0 seconds. 

Real Correlation time (seconds) for changes 

in the cross -track gravitational 
accelerations. Default =1.0 seconds. 

Real Correlation time (seconds) for changes 

in the along-track gravitational 
accelerations. Default =1.0 seconds. 
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3.3.16 MBIAS 



The MBIAS card is used to initialize the standard deviations for station 
measurement biases. If omitted, no biases will be adjusted. 


MBIAS 101 0 0,001 0,000001 
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Variable 


Type 

Description 

1 

Alphanumeric (A8) 

Keyword MBIAS. 

2 

Integer 


Station number (if 0, this card will 
apply to all stations). 

3 

Integer 


Number of biases (1 or 2). Default=l. 

4 

Real 


k priori standard deviation for the 
first measurement bias (kilometers for 
range). No default. 

5 

Real 


A priori standard deviation for the 
second measurement bias (km/sec for 
range rate). No default. 

Note: 

Note (1) on 

the CSTA 

card also applies. See also note on MSIG card. 
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3.3.17 REFRAC 






..tij r 


The REFRAC card is used to initialise the standard deviations for station 
refraction parameters. If omitted, no refraction parameters will be adjusted. 

/REFRAC 101 0 0, 0005 ,,,, 
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MCCC-90«I 


Variable 

Type 

Description 

1 

Alphanumeric (A8) 

Keyword REFRAC. 

2 

Integer 

Station number (if 0, this card will 
apply to all stations). 

3 

Integer 

If this number is greater than 0 on 
any REFRAC card, then the program 
will use the refraction correction 
on the measurement tape as the 
multiplier of the refraction parameter. 
Default = 0. 

4 

Real 

k priori standard deviation for the 
measurement refraction parameter (see 
Note 2 ) . No default. 
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Note; (1) Note (T) on the GSTA card also applies. 

(2) If variable 3 Is zero for all REFRAC cards, then the 
residual refraction error Is assumed to be of the form: 


K/sinE (range) 

'<<-5?^ - '""Se «te) 


where K 1s the refraction parameter and E is the 
elevation of the ray. Thus the a pHoH sigma on K has 
units of kilometers. If variable 3 is greater than zero 
for any REFRAC card, then K will be treated as the multi- 
plier of the total refraction correction which is available 
on the measurement tape for each measurement. Thus K is 
dimensionless. 
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3.3.18 STATEl 


ORIGINAL PAGE IS 
OF POOR QUALITY 


The STATEl card is 
the first three orbital 


used to initialize the standard deviations for 
states. 
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Variable Type 


Description 


1 

2 


3 

4 

5 

6 


Alphanumeric (A8) Keyword STATEl. 

Integer Defines type of orbital elements: 1 = 

cartesian elements; 2 = Keplerian elements; 
3 = spherical elements. Note, this must 
match type on ORBIT file. Default = 1 

Integer Not used. 

Real A priori standard deviation for satellite 

State 1 (km). Default = .01^^^. 

Real A priori standard deviation for satellite 

state 2. Default = .01^^^ 

Real A priori standard deviation for satellite 

state 3. Default = .01^^^. 


Note: (1) Units are kilometers, radians or dimensionless. 
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STATEIQ 


3.3.19 STATEIQ 

The STATEIQ card is used to input the state noise on the satellite 
position in earth-centered inertial cartesian coordinates (the program 
makes the transformation to other types of orbital elements). 
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Variable Type 


Descri pti oti 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword STATEIQ. 

Not used. 

Not used. 

Square root of state noise spectral 
density for the x component of satel- 
lite position (in km/sec^^^). Default = 0. 
Square root of state noise spectral 
density for the y component of satel- 
lite position (in km/sec^^^). Default = 0. 
Square root of state noise spectral 
density for the z component of satel- 
lite position (in km/sec ' ). Default = 0. 
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STATE2 


ORIGINAL PAGE \S 

3.3.20 STATE2 OF POOR QUALITY 

The STATE2 card is used to initialize the r.yf.lard deviations for 
orbital states four, five and six. 
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Variable 

. _ Type. 

Description 

1 

Alphanumeric (A8) 

Keyword STATE2. 

2 

Integer 

Not used. 

3 

Integer 

Not used. 

4 

Real 

A priori standard deviation for satellite 
state four. Default = .0001^^^. 

5 

Real 

A priori standard deviation for satellite 
state five. Default = .0001^^ \ 

6 

Real 

A priori Standard deviation for satellite 
state six. Default = .0001^^^. 


Note: (1) Units are km/sec or radians. 
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STATE2Q 


I 


3.3.21 STATE2Q 


ommi m 

OF. POOR QUALITY 


The STATE2Q card is used to input the state noise on the satellite 
velocity in earth-centered inertial cartesian coordinates (the program 
makes the transformation to other types of orbital elements). 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword STATE2Q. 

Not used. 

Not used. 

Square root of state noise spectral 
density for the x component of satel- 
lite velocity (in km/sec^'^^). Default = 0. 
Square root of state noise spectral 
density for the y component of satel- 
lite velocity (in km/sec^^^). Default = 0. 
Square root of state noise spectral 
density for the z component of satel- 
lite velocity (in km/sec^'^^). Default = 0 
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THRUSTl 


ORIGINAL' PAGE fS 

3.3.22 THRUSTl OF POOR QUALITY 

The THRUSTl card is used to initialize the standard deviation for 
the perturbing components of the first thrust segment. If omitted, the 
perturbing thrusts will not be adjusted. 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword THRUSTl . 

Not used. 

Not used. 

4 priori Standard de'. ■nation for the 
vertical component of thrust acceleration 
(km/ sec ). No default. 

4 priori Standard deviation for the cross- 
track component of thrust acceleration 
(km/sec ). No default. 

4 priori Standard deviation for the along 
track component of thrust acceleration 
(km/sec^). No default. 
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THRUSTIQ 


3.3.23 THRUSTIQ 


ORIGINAL' PAGE fS 
OE POOR QUALITY 


The THRUSTIQ card is used to input the state noise for the perturbing 
components of the fi rst thrust segment. 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword THRUSTIQ. 

Not used. 

Not used. 

Square root of state noise spectral 
density for the vertical component of 
thrust acceleration (km/sec^'^^) . 
Default = 0. 

Square root of state noise spectral 
density for the cross-track component 
of thrust acceleration (km/sec^'^^) . 
Default = 0. 

Square root of state noise spectral 
density for the along-track component 
of thrust acceleration (km/sec^^^). 
Default =0. 
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3.3.24 THRUST2 


PmE IS 
OF POOR QUALITY 


The THRUST2 card is used to initialize the standard deviations for 
the perturbing components of the second thrust segment. If omitted, 
the perturbing thrusts will not be adjusted. 
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Variable 

1 

2 

3 

4 


Type 


Alphanumeric (AS) 
Integer 


Description 


Integer 

Real 


Real 


Real 


Keyv/ord THRUST'c. 

Not used. 

Not used. 

A priori, standard deviation for the 
vertical component of thrust acceleration 
(km/ sec ). No default. 

A prioH standard deviation for the cross- 
track component of thrust acceleration 
(km/sec^). No default. 

A priori standard deviation for the along 
track component of thrust acceleration 


(km/sec^). No default. 


I 
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3.3.25 THRUST2Q 


ORIGIWAL’ PAGE !S 
PE POOR QUALITY 


The THRUST2Q card is used to input the state noise for the perturbing 
components of the second thrust segment. 
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Variable Type 


Description 


1 

2 

3 

4 


5 


6 


Alphanumeric (A8) 

Integer 

Integer 

Real 


Real 


Real 


Keyword THRUST2Q. 

Not used. 

Not used. 

Square root of state noise spectral 
density for the vertical component of 
thrust acceleration (km/sec^^^). Default = 0. 
Square root of state noise spectral 
density for the cross- track component of 
thrust acceleration (km/sec^^^). Default = 0. 
Square root of state noise spectral 
density for the along-track component of 
thrust acceleration (km/sec^^^). Default = 0. 
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4.0 PROGRAM OUTPUT 


The printed output from PREFER consists of the filter estimates at 
the end of each mini -batch, a summary of the measurement residual statis- 
tics for each mini-batch, the smoother estimates at the corresponding times, 
and optional printer plots of the filter and smoother position standard 
deviations. PREFER also produces an ORBl File (smoother output) on unit 
70 and an optional file (unit 71) containing the smoother covariance matrix 
in ECI coordinates. 

The filter summary includes output in both earth-centered inertial 
or local coordinates at the current time and epoch elements (cartesian, 

Keplerian or spherical). The ECI/HCL output consists of the nominal 
trajectory, the estimated correction and the estimated total. The fil- 
ter output for the epoch elements is the estimated correction to the 
nominal elanents. Also printed are the estimates for the remaining state 
parameters and the computed standard deviation for the estimate. 

The measurement residual statistics for each mini-batch are printed 
before the corresponding filter summary. Two lines are printed for measure- 
ments which occur in pairs (e.g., GPSPAC, or ground range and range rate). 

The first line for each GPS satellite lists the statistics for the pseudo 

range measurements while the second line corresponds to the pseudo delta- 
range measurements. Ground tracking can be easily distinguished from satellite- 
to-satellite tracking because GPS satellite numbers are printed as negative. 
Included in the residual statistics is the "WEIGHTED SIGMA". This is computed 

as J (z(z?/o 2 ))/n where a is the residual sigma ( 02 - = / HPH^+R ) 
i i ^i 

as computed by the filter. If all the filter models are correct, the expected 
value of the WEIGHTED SIGMA should be approximately equal to 1 and the 
deviation from 1 should be approximately l//2rT . 

Between the filter sumnaries may be other messages which indicate 
when pass parameters are added or removed. Also, any measurement residuals 
which exceed the editing or printing thresholds will be printed. 
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A summary of the measurement residual statistics is printed at the 
end of the filter printout. The following five quantities are printed: 

(a) Total number of measurements: n 


(b) Weighted sigma: 

n 0 0 

(c) i zVoi 

i-1 ’ ^i 

(d) E m{al ) 
i»l ^i 

(e) -i Z (z?/a^ al ) . 

The sum of weighted residuals (quantity c) should be chi-square 
distributed with a mean of n and a standard deviation of . If the 
value printed deviates greatly from n ± , then modeling errors are 

suspected. Quantity d is the sum of the natural logarithm of the 
analytically computed residual variances. Quantity e is the log likeli- 
hood function (plus a bias) and is computed as the negative of one-half 
the sum of quantities c and d . This is the quantity which would be 
maximized in a maximum likelihood estimation program [3]. Thus, the log 
likelihood can be used as a metric to compare runs made with different 
model assumptions. Notice that this procedure attempts to minimize the 
sum of weighted residuals and the residual variances: it is not acceptable 
to minimize just the sum of weighted residuals since this can be done by 
making a_ large. 


/ (I(2?/4 ))/" 

i ’ ^i 
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The smoother output is printed backward in time after the filtering is 
completed. The smoother summary is very similar to the filter summary and 
also Includes the filter estimate at each point. 

2 

All output is in units of kilometers, km/sec, km/ sec , radians, 
milliseconds. Hz/megahertz or dimensionless. 

4.1 Format of ORBl File (Unit 70) 

The smoothed ECI positions and velocities are written on Unit 70 in 
the format of a GTDS ORBl File. This will be generated in the same 
coordinate system (TOD or 1950) as the input ORBIT file. The starred 
quantities are not actually written on the file. 


Mcono 

BKomo 

size 

iVTlS 

NAME, csscRiPTiCN, rrc. 

1 

2bOO bycas 

1-8 

TlTLE(l ) : 

ORBl Idsnelflsr ■ 76796291. 

08B1 

txMder 


!J-16 

lOSAX: 

SaCsUlca ID nuafasr 



L7-24 

tloe used 




25-32 

STOKBIU) : 

Start of sphsmsris (.YYMMDO.) 



33-40 

TITLE<5) ; 

Day count of year for start of 
sphsmsris 



41-48 

TlTLii(6); 

Sseonda of day for start of 
sphsmsris 



49^56 

ETOBBIU): 

Exxd tima of sphsmsris (THIMDO.) 



57-84 

T17LE(b) : 

Day count of year for end of . 

sphsmsris 



65-72 

HTLE(9) : 

Seconds of day for end of 
sphsmsris 



73-80 

TITLE(iO) ; 

Tims interval between ephemaris 
points (sac.) 



81-152 

HEADER CI,1) 

; 1-1,9 





Title of current run 



153-160 

ICS2IT: 

Central body indicator 
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r 


i 


nftCQHD 

MGOnO 

SlZi 

•YTIS 

NAMC, oeSCniPTION, ETC. 

1 

ORBl 

header 

2800 bytes 

161-208 

209-216 

Spares 

YMDR£F: 

Year» monthy and day of referencec 
elsM of true system 



217-224 

TITLZ(28) : 

Day count of year for reference 

data 



225-232 

GHA: 

Greenwich hour angle at epoch 



233-632 

Spares 




633-MO 

641-648 

649-656 

CSUBDZ: 

AREA: 

SCMASS: 

Drag coefficient 

Area of spacecraft In cm3 

ttosa of spacecraft in grams 



657-800 

Spares 




801-808 

TITLEdOlJ ; 

Tima of epoch from reference date 
in seconds /806. 81242 



809-816* 

TITLEU02) : 

Ratio of semimajor axis Co mean 
equatorial radius 



817-824* 

AEINTU) : 

Eccentricity 



825-832* 

0BLINT(5): 

True anomaly at t^: rad. 


i 

1 

833-880 

pvmr(i) : 

1-1,6 

Carc^ian elements for incegratio: 
initialization 


1 

1 

88l-aS8 

889-896 

SPINT(5) : 
®PUITC6) : 

Magnitude of radius vector 
(km) 

Magnitude of velocity vector 
(km/sec) 

Sum of mean anomaly and argument 
of perlfocus 


j 

897-904* 

T1T1E<U3) : 



905-912* 

AEINTC6) : 

Maan anomaly 



913-920 

Spare 




921-928* 

lAEINICS): 

Augment of perlfocus 



929^936* 

AEINT(3) : 

Inclination 



937-944* 

A£IST(4) : 

Longitude of ascending node 


i 









nicono 


1 

ORBl 

header 


nccoRO 

SIZf 



2800 bytes 945-952* 


95>960* 


NAMI, OISCRIFTION, 6TC. 


TITLE (119): Vertical flight path angle alnua 

r/2 

0BLIMT(4): Mean motion at t^: rad. /DOT sec 


981-968* OBLIMT(l): Eccentric anomaly at t^: rad. 

969-976* 0BLmT(8): Time rate of change of argument of 

perigee at C^: rad. /DOT sec. 

977-984* 0BL1NT(9): Time rate of change of right ascen- 
sion of ascending node at t.: rad. 

DOT sec. ^ 

985-992* 0BLHIT(2): Period at DOT sec. 


993r-1000 


1001-1008 


1009-1520 


IBLIMT(6): Perigee height c^: km 

iBLIllt(7): Apogee height at km 


1521-1528 

1529-1536 


1537-1544 


1545-1552 


1553-1560 


1561-1568 


Tear of elements epoch 
Month of elements epoch 
Day of elements epoch 
Hour of elements epoch 
Minute of elements epoch 
Seconds of elements epoch 


1569-1592 


1593-1600 


1601-2800 


[TIE (200): Orbit theory Indicator 

- 2, Cowell or time-regularized 
Cowell 

■ 3, all odiers 
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nKono 

RKCOnO 

SIZE 

•ms 

namI, DiscmmoN, era 

All cxctpe 
first and 
last two 

2800 bytas 

1-8 

9-16 

DATA(l) : 
DATA(2 ) i 

Data, of first aphsmaris point 
(TEHMDD.) 

Day count of yaar for first 
aphemaris point 



17-24 

DATA(3) : 

Saconds of day for first aphamaris 
point 

ORBl data 


23-32 

DATA(4) : 

Tiaa Interval batwaan data points 
(sac.) 



33-40 

DATA(5) ; 

Day count of year for epoch 



41-86 

PVINT(I) : 

I • 1,6; first position vector 
(kn) and velocity vector (kn/sec) 



87-2440 

PVINT(I) ; 

I » 1,6; J = 2,50 
Position and velocity vector sets 
for data points 2-50. For less 
than 50 points, the first invalid 
set is .9999999999999999 x lO^® 
(units ka and kn/sec) 


1 

2441-2800 

Spares 


Last two 
ORBl and 
santlnols 

2800 bytas 

1-2800 

ESNINE(I) 

I = 1,350 

End of file sentinels, with the 
value 0.9999999999999999 x lOl® 
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4.2 Format of Smoother ECI Covariance File (Unit 71) OF POOR QUALiyy 

If requested, the entire smoother ’covariance matrix is written on 
unit 71 where the entries corresponding to the orbital elements are 
transformed from epoch osculating elements to current ECI cartesian 
elements. This file consists of two types of records which are repeated 
until the end of the trajectory is reached. Notice that the first 
record (in a pair) contains a variable which defines the length of the 
second record. 

Record Type 1 


Word 

No. 

1-2 

3 

4 

5 

6 


Type 

Double Precision 

Integer 

Integer 

Integer 

Integer 


Description 

Time from epoch (sec). 

Number of adjusted states (n) at current time. 
Length of covariance matrix 

Numeric label of 1st state. 

Numeric label of 2nd state. 


104 Integer 


Numeric label of 100th state. 
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Record Type 2 


This contains the ECI smoother covariance matrix stored as upper 
triangular by columns. 


Word 

No. 

1-2 

3-4 


Type 

Double Precision 
Double Precision 


Description 

Variance of 1st state. 

Covariance between 1st and 2nd states. 


(n» (n+l)-l) Double Precision Variance of nth state, 
-{n* (n+1)) 


Notes; (1-) Because of numerical problems, it is possible for the state 
variances to become negative. 

(2) The numeric labels of the states are defined as: 

Label Description 

1-6 state orbital elements (ECI) 

7 d?ag coefficient 

8 perturbing gravitational accelerations (HCL) 

9 acceleration of first thrust segment (HCL) 

10 acceleration of second thrust segment (HCL) 

11 host clock timing error (msec) 

12 host clock drift rate (msec/sec) 

13 altimeter bias 

100000 station bias 

200000 station refraction parameter 

300000 station position error (ECI) 

400000 GPS position error (HCL) 

500000 GPS timing error (msec) 

The last 5 digits in each number are reserved for the 
station or satellite number. 
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5.0 EXAMPLE (LANDSATl 


Orbit - 700 km altitude, e = .0001, 98.3® inclination, 720 minutes 
(7.3 revolutions) 

Simulated Tracking Data - 5 ground stations measuring range and 

integrated range rate. See Table 5.1 for error parameters. 

Nominal Trajectory - Obtained by using GTOS to least squares fit the 
simulated tracking data over 12 hours. See Table 5.1 for 
differences between models used in simulating data and models 
used in batch orbit determination. 

Table 5.1 Error Parameters of Example 



Parameter 

Model 

Description 


Simulated Data 

Batch Orbit Determination 

data noise (range) 

1 meter 


1 meter (weighting) 


data noise (range rate) 

0.1 cm/sec 


0.1 cm/sec (weighting) 


gravitational field 

GEM9 (20,20) 


WGS72 (16,16) 


measurement bias (R) 

5 meters 


— 


measurement bias (ft) 

0.1 cm/sec 


— 


station position errors 

±5m each component 




drag coefficient 

2.0 


1.5 


density model 

Flux table #150 


Flux table #150 


solar radiation pressure 

I 

yes 

no 

ref racti on-tropospheri c 

yes 

yes 

ref racti on-i onospheric 

yes 

yes 

polar motion modeled 

no 

yes 











RSS POSITION ERROR (METERS) 
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A priori standard deviations used in PREFER 


x,y,2 

• • • 

x,y,2 

gravitational acceleration 

station measurement biases 
station refraction 
station position errors 


- 20 meters 


2 cm/sec 

LO 
.2 


2.4xlO’^(H), 0.7xl0"®(C), 1.0xl0"^(L) 


m/sec 

5 meter (range), 0.1 cm/sec (range rate) 
50 cm @ 2 enith 
5 meters (each component) 


State noise spectral density - none on orbital elements 


Markov process standard deviations 

gravitational accelerations - 2.4x10"® (H), 0.7xl0"®(C), 1.0xi0"®(L) 

2 

m/sec 


Time constants of Markov processes 

gravitational acceleration - 1200 (H), 200 (C), 200 (L) seconds 

Results - See Figure 5.1 and Table 5.2 


Table 5.2 Position Errors for Smoothed Trajectory 


component 

RMS 

maximum 

radial 

3.8 

8.1 

crosstrack 

7.0 

13.1 

alongtrack 

11.4 

24.9 

total 

13.9 

26.0 


The nominal trajectory has large errors (e.g. 89 meters) in the data 
gaps hear the end of the orbit. This occurs when only one station (#31) is 
tracking. In the periods when two or more stations are tracking, the error 
in the nominal trajectory is less than 32 meters. 
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The smoothed output of PREFER has errors which are less than 21 
meters during the period when two or more stations are tracking. In the 
last data gap at the end of the trajectory, the peak error is 26 meters 
(mostly alongtrack). Notice that the total error in the PREFER output 
tends to be relatively smooth while the error in the nominal trajectory 
(GTDS output) fluctuates wildly; there are 4 short periods of time when 
the nominal trajectory is more accurate than the PREFER output but the 
PREFER output is generally much superior. 

These results for PREFER were obtained from one of the best runs in 
a series of runs in which a priori variances, state noise variances, time 
constants, and the selection of adjusted parameters were varied. Several 
other runs yielded likelihood functions and ephemeris errors similar to 
those of the run listed. It was found that the results are not very 
sensitive to the exact form of the gravitational acceleration model j_f 
the parameters of that model are reasonable. Maximization of the likeli- 
hood function can generally be relied upon to select the run with the 
vptimum model . 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

THE PRECISION RECURSIVE ESTIMATOR FOR EPHEMERIS REFINEMENT IPREFER) VERSION (IBM) 82.001 

OUTPUT FROM ORBIT FILE FOR SATELLITE LS-D 
EPOCH OATE/TIME - 811001. 195500.0000 

START OATE/TIME - 811001. 195500.0000 END DA1E/TIME = 811002. 200000 OOOO 

SATELLITE AREA - 20.000 M-*2. MASS - 1700.00 KG, CD - 1.500 

ORBIT FILE MAS CREATED USING 1950 COORDINATES 

EPOCH KEPLERIAN EteMENTS • 7078.032 0.000099 I.71A872 5.846084 0 520692 2.620901 

EPOCH CARTESIAN ELEMENTS • -6413.078 2996.587 -0.687 0.45640 0.97591 -7.42596 


oo 
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LISTING OF CARO INPUT 


'■wm 


f 


PRINT 

2203 

0 

.O 

,o 

.0 

o 

5 

EARTH 

0 

0 

6378.1660 

298 .30000 

398600.64 

.299790406 

c 

STEP 

O 

o 

120.00000 

43200.000 

120.00000 

0 


ELCUT 

o 

o 

2.0000000 

.O 

.O 

.o 

STATE 1 

0 

o 

.200000000-01 

.200OO0QOD'OI 

.200000000-01 

O 

Q 

STATE2 

0 

o 

.200000000-04 

.200000000-04 

.200000000-04 

O 

1 

GRAV 

0 

0 

. 240000000-08 

.700000000-09 

. 100000000-OB 

.0 

GRAVO 

o 

o 

.240000000-08 

.700000000-09 

.100000000-08 

O 


GRAVT 

0 

0 

1200.0000 

200.00000 

200.00000 

.0 


MBIAS 

0 

2 

.5O0G0O00O-O2 

. lOCOOOOOO-OS 

.0 

0 

Cl 

REFRAC 

o 

o 

.500000000-03 

.o 

.0 

0 

2 

CSTA 

0 

o 

.500000000-02 

.500000000-02 

.500000000-02 

0 

MSIG 

4 

2 

. 100000000-02 

. 100000000-05 

.O 

O 


MSIG 

13 

2 

. 100000000-02 

. 1000000CD-05 

.O 

0 


HSIG 

30 

2 

. 100000000-02 

. 100000000-05 

.O 

O 

Xa 

MSIG 

31 

2 

. 100000000-02 

.100000000-05 

.O 

■O 

a 

MSiG 

32 

2 

. 100000000-02 

.100000000-05 

.0 

.0 

1 

i 




1 
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SUMMARY OF STATE PARAMETERS 



LABEL 

LABEL 

UNITS 

A PRIORI SIGMA 

SQRTI SPECTRAL OEN. ) 

TIME CONSTANT 

NOISE SIGMA 

1 

1 

X 

KM. 

.20000-01 

.0 



2 

2 

Y 

KM. 

.20000-01 

.0 



3 

3 

2 

KM. 

.20000-01 

.0 



4 

4 

XDOT 

KM/SEC 

.20000-04 

0 



5 

5 

YDOT 

KM/SEC 

.20000-04 

.0 



6 

6 

ZDOT 

KM/SEC 

.20000-04 

O 



7 

8 

GRAV 

KM/S«*2 

.24000-08 

.97980-10 

1200. 

24000-08 

6 

8 

GRAV 

KM/S**2 

.70000-09 

.70000-10 

200.0 

70000-09 

9 

8 

GRAV 

KM/S*»a 

. 10000-08 

. 10000-09 

200 O 

, 10000-08 

10 

100000 

MBIAS 1 

KM. 

.50000-02 

.0 



It 

200000 

MBIAS 2 

KM/S 

10000-05 

.0 



12 

300000 

REFRACT 

KM.,- 

.50000-03 

.0 



13 

400000 

S7A POS 

KM. 

.50000-02 

.0 



14 

400000 

STA PQS 

KM. 

.50000-02 

.O 



15 

400000 

5TA POS 

KM. 

.50000-02 

.0 




NOTE: FOR A FIRST ORDER MARKOV PROCESS. OUTPUT SIGMA - SORTlENPUT SPECTRA!. DENSITVTIME CONSTANT/3) 


O 


o o 

m 2 

■o S 
o z 
o ^ 

70 n 

<o 
c > 
■> o 
r; ns 

-< c: 
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SUMMARY QF GROUND Sr^tTIONS ( 5) 

station posit IONTEFCI input HeASUREMENT SIGMAS 



STATION NO. 

X 

Y 

z. 

SIGMAI 

S1GMA2 

I 

4 

4847.8243 

-353.3178 

4117. 1445 

. 10000-03 

. 10000-05 

2 

31 

-2282.4950 

-1453.3850 

5756.7222 

. 10000-02 

.10000-05 

3 

30 

-4447. 490G 

2676.8624 

-3695.2795 

. 10000-02 

, 10000-05 

4 

13 

1129.8039 

-4833. 1549 

3992.2086 

. 10000-02 

.10000-05 

5 

32 

1263.4446 

-6255 0397 

-68.8004 

. 10000-02 

. 10000-05 










HOTES ON RESIDUAL STATISTICS 

t| If ALL MODELS ARE CORRECT. THE MEIGHTEO RESIDUAL RMS SIttULD EQUAL I. UITM A SIGMA OF 1,/S0RTC2*Nl 
2) UNITS OF RESIDUALS ARE KM (ALL TYPES EXCEPT DOPPLER) OR DIMENSIONLESS (DOPPLER) 


APPROXIMATE REGION SIZE » 3B32 16 BYTES (ASSUMING TWO 7294-BVTE BUFFERS ARE USED FOR UNITS 70 AND 7t) 

REGION SPECIFIED FOR THIS UOB * 405S04 ARRAYS REQUIRE 26360 BYTES 

NOTE: PREFER MILL USE ALL THE CORE ALLOCATED FOR THE JOB. IGNORE REGION PRINTED IN SYSTEM MESSAGES 
5288 4496 6040 5288 5288 79689 78351 76429 7508 1 73545 


FILTER OUTPUT AT TIME 120.00 


NOMINAL TRAJECTORY (ECI)- - 

6306.5691 

3089. 1352 

-889.3853 

1,31634 

0.56448 

-7.36562 




ESTIM. TOTAL (ECI) > -6306.5691 

3089. 1352 

-889,3853 

1,31634 

0,56448 

-7 36562 




ESTIM. CORRECTION (HCL) - 

0.0 

0.0 

0.0 

0.0 

0-0 

0.0 




ESTIM. SIGMA (HCL) 

0.02045 

0.01998 

O 02001 

0.00002 

O 00002 

O 00002 




PARAMETER LABELS * X 

V 

Z 

XDOT VOOT ZOOT <JRAV 

GRAY 

GRAY 


STATION/SATELLITE 

O 

o 

O 

O 

O 

0 

O 

O 

O 

FILTER STATE (EPOCH) - 0.0 

0.0 

0.0 

0.0 

O O 

0.0 

0.0 

0.0 

0 O 



FILTER SIGMAS (EPOCH)- 0.02000001 0.02000000 0.02000000 0,00002000 0.00002000 0.00002000 0.00000000 0,00000000 0.00000000 


FILTER OUTPUT AT TIME 240.00 

NOMINAL TRAJECTORY (ECl)- -6098.0154 3131.6999 -1763.6533 3.15483 

ESTIM. TOTAL (ECI) - -6098.0154 3131.6999 -1763.6533 2.15483 

ESTIM. CORRECTION (HCL) « 0.0 0.0 O O 0.0 

ESTIM. SIGMA (HCL) • 0.02161 0.01993 0.02021 0.00002 


O 14399 -7 18580 

0.14399 -T. 18580 

0.0 O O 

O 00002 O 00002 


to 

IN3 


PARAMETER LABELS 

* X 

Y 

z 

XDOT 

YOOT 

ZOOT 

GRAY 

GRAY 

STATION/SATELLITE 

« 

o 

o 

O 

O 

O 

O 

O 

FILTER STATE (EPOCH) 

- O 0 

0.0 

0.0 

0.0 

O O 

0.0 

0.0 

0.0 


FILTER SIGMAS (EPOCH)- 0.02000009 0.02000003 0.02000001 0.00002001 0.00007000 O 00002000 O OOOOOOOO 0.00000000 O.OOOOOOTO 


FILTER OUTPUT AT TIME 360.00 


NOMINAL 

TRAJECTORY (ECI)- 

-5790.8142 

3123.6042 

-2609,3150 

2.95823 

-0.27871 

-6 88950 

ESTIM. 

TOTAL (ECI) 

-5790.8142 

3123 6042 

-2609.3150 

2.95833 

-O 37871 

•6 88950 

ESTIM 

CORRECTION (HCL) - 

0 O 

0.0 

0.0 

0.0 

O O 

0 O 

ESTIM. 

SIGMA (HCL) 

0.02313 

0 01985 

0 02107 

0,00003 

O 00002 

O 00002 


parameter labels - X 

STATION/SATELLITE 

Y 

O 

2 

O 

XDOT 

O 

VDOT 

O 

O 

FILTER state (EPOCH) » 0,0 

0-0 

0.0 

0.0 

0.0 



FILTER SIGMAS (EPOCH)* 0.02000044 O 02000013 O 02000(^7 0.00002001 0.000020CX> 


GRAY 


O O 


0.0 


0.0 


0.0 


0.00002000 0.00000000 0.00000000 0.00000000 


o o 

“rj 

Tj 1 
O S 
0 5 » 
w n 

(O T3 

o 

ET 


FILTER OUTPUT AT TIME 480.00 

NOMINAL TRAJECTORY (ECI)« -5389.9667 3064,9966 -3412.6691 

ESTIM. TOTAL (ECI) * -5389.9667 3C64.9966 -3412.6691 

ESTIM. CORRECTION (HCL) * 0.0 0.0 0.0 

ESTIM. SIGMA (HCL) * O.0245S 0.01975 0.02319 


3.7134S 

3.71349 

0.0 

0.00003 


-O 69674 
-O 69b/4 
0.0 

O 00002 


-6 48164 
-6.48164 

0.0 

O 00003 


PARAMETER LABELS 

» X 

Y 

2 

XDOT 

YOOT 

200T 

GRAY 

GRAY 

STATION/SATELLITE 


O 

O 

O 

O 

O 

O 

O 

FILTER STATE (EPOCH) 

* 0.0 

0.0 

0.0 

0.0 

0.0 

O O 

0.0 

0.0 


FILTER SIGMAS (EPOCH)* 0.02000135 0.02000041 0.02000027 0,00002003 0.tXXK)200l 0.00002000 0.00000000 0,00000000 O-OOCOOOOO 
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FILTER StCSMAS fEPOCH|> 0.0202C99R 0.02022S77 0.02229970 O 00002219 O 00002064 O 00002039 0.00000000 0.00000000 0.<X)000000 




FILTER OUTPUT AT TIME 28BO.OO 


NOMINAL TRAJECTORY 

(ECD- 

6421 1216 

-2899. 1272 

-601 8361 

0. I30S5 

-1 24688 

ESTIM. TOTAL (ECI) 


6421.1218 

-3839. 1272 

-601.8361 

0.13095 

-1 24688 

CSTIM. CORRECTION 

(HCL) - 

0.0 

0.0 

0,0 

0-0 

0.0 

ESTIM. SIGMA (HCL) 

» 

0.07318 

0.03457 

0.36930 

0.00018 

O 00018 


7 4064S 
7.40645 
0.0 

O 00013 


PARAMETER LABELS 
STAT lON/SATELLl TE 
FILTER STATE (EPOCH) 


XDOT 


YDOT 


0.0 


0,0 


0.0 


0,0 


0.0 


0.0 


0.0 


0,0 


0.0 


FILTER SIGMAS (EPOCH)- 0.02029040 O 02025363 0. 03371159 0,00002349 0.00002073 0.G0Ci02040 0.00000000 0.00000000 0,00000000 


filter output at time 3000.00 


NOMINAL TRAJECTORY (ECI)- 

6384.6691 -3024.8148 

383.4234 

-0.73767 

-0.84507 

7.42769 

ESTIH. 

TOTAL (ECI) 

6394.6691 -3024.8148 

389.4234 

-0.73767 

-O 94S07 

7 42769 

ESTIM. 

CORRECTION (HCL) - 

0.0 0.0 

0.0 

0.0 

O O 

O O 

ESTIM. 

SIGMA (HCL) 

0.12520 0.03706 

0.26953 

0.00020 

0.00021 

0 00007 


PARAMETER LABELS 

X 

Y 

2 

XOOT YOOT 

ZDOT 


GRAV 

GRAV 

GRAV 

STATION/SATELLITE 


O 

o 

O O 

O 

O 

O 

0 

o 

FILTER STATE (EPOCH) - 

0.0 

0.0 

0.0 

0.0 0.0 

O O 


0.0 

0.0 

0.0 

FILTER SIGMAS (EPOCH)- 

0.02028880 0.02027626 0.02316704 0.00002280 0.00002081 0.00iX>204l 

0.00000000 

0.00000000 

0.00000000 


FILTER output AT TIME 3130.00 


NOMINAL TRAJECTORY (ECI)« 

6244.5604 -3101.3954 

1175 9713 

-1.59439 

-0 ^<2955 

7 33S03 



ESTIM. TOTAL (ECI) 

6344.5604 -3101.3954 

1175.9713 

-1.59439 

-O 63955 

7 33803 



ESTIM. CORRECTION (KCL) - 

0.0 0.0 

o.d 

0.0 

O G 

o e 



ESTIM. SIGMA (HCL) 

0.17785 0.07387 

0.35398 

0.00020 

O ^24 

0 oocaoa 



PARAMETER LABELS • X 

Y 2 

XDOT YDOT 20QT GRAV 

GRAV 

GRAV 

STATION/SATELLITE 

O O 

O 

O 

0 


o 

O O 

FILTER STATE (EPOCH) • 0.0 

0.0 0.0 

0.0 

0.0 

0.C 

0.0 

0.0 

0.0 


FILTER SIGMAS (EPOCH)- 0.02029544 0.03029941 0.03365193 0.00003313 O.O0OD209O 0.00003043 O OOOOOOOO 0^00000000 0.00000000 


o o 

Jr «==*» 

o ^ 
c ^ 
» r % 


c: 


< 


m 


FILTER OUTPUT AT TIME 3340.00 


NOMINAL TRAJECTORY 

(Ecn- 

6003.0835 

-3127.6322 

2043.3798 

-3.42483 

-O 00716 

7. 10915 

ESTIM. TOTAL (ECI) 

» 

6003 0835 

-3137.6332 

3043.3798 

-2 42483 

-0 00716 

7 10915 

ESTIM. CORRECTION 

(HCL) - 

O O 

O O 

O O 

O O 

O O 

O O 

ESTIM. SIGMA (HCL) 

s 

0.23597 

0. 11701 

0.21887 

0.00018 

0 00026 

O OOOlO 


PARAMETER LABELS 

» X 

Y 

2 

XOOT 

YDOT 

200T 

GRAV 

GRAV 

GRAV 

STATION/SATELLITE 

» 

O 

O 

O 

O 

O 

O 

O 

O O 

FILTER STATE (EPOCH) 

» 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 


FILTER SIGMAS (EPOCH)- 0.02030071 0.03032182 0.02416100 0.00003347 0.00002099 O 00002043 0.00000000 0.00000000 O OOOOOOOO 


NEW PASS PARAMETERS ADDED AT TIME 
%%%% RESIDUAL STATISTICS FOR TIME 


3420.03. STAT30N/SAT 4. N » 15 

3300.00 TO 3420.00 








FILTER OUTPUT AT TIME 3360.00 


BusiSEssAND TeamwGicAL Systems, Inc 



V / 



NOMINAL TRAJECTORY (EC1>« 

5664.1891 

-3103.1151 

2877.5419 

-3.21570 

O 41519 

6 77470 

estim. 

TOTAL (ECI) 

5664 . I860 

-3103. 1131 

3877.5292 

-3,21570 

0.4IS19 

6 77471 

ESTIM. 

CORRECTION (HCL) * 

0.00379 

0.00244 

-0.01243 

-0.00000 

-0.00000 

O 00001 

ESTIM. 

SIGMA (HCL) 

0.00989 

0.01951 

0.01065 

0,0000% 

O 00002 

O 00002 


parameter labels ■ X V 2 XOOT VOOT 200T CRAV GRAV GRAV MB I AS I 

STATION/SATELLXTE - O O O O O O O O 0 4 

FILTER STATE (EPOCH) » 0.00053431 O.OOtSO I IR-0. 00377666-0. 00000354 0.00000254 0.0000017 1-0. 0000CMX)d-O.5bo606bO 0. 00000000-0. OOOOatsi 
FILTER SIGMAS (EPOCH)* 0.01611238 0.01901381 0.02153840 0.00002093 0.00001928 0.00001360 0.00000000 0.00000000 0.00000000 0.00499854 


parameter labels 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


- MBIAS 2 REFRACT STA POS STA POS STA POS 
• 4 4 4 4 4 

•-0.00000000 0.00002852-0.00001858 0.00002447 0.00014427 
» 0.00000100 0.00049182 0.00499965 0.00499940 0.00497017 


ftft RESIDUAL STATISTICS FOR TIME 
STATION 4 TYPE tO N* 

STATION 4 TYPE 13 N« 


3430.00 TO 3540.00 

12. MEAN* .74960-07 SIGMA* .23060-05 WEIGHTED SIGMA* 

Tir. MEAN- -.17700-03 SIGMA* .78460-03 WEIGHTED SIGMA* 


1.044 

.6938 


FILTER OUTPUT AT TIME 3480.00 


NOMINAL TRAJECTORY (ECU* 

5233.4213 -3(^73 3663 

3664 9082 

-3.95398 

0 830S6 

6 33031 

ESTIM. 

TOTAL (ECI) 

5233.4172 ^3028.2577 

3664.8985 

-3.95398 

0.83056 

6.33031 

ESTIM. 

CORRECTION (HCL) « 

0.00160 

-0.r00063 

-0.01052 

-0.00001 

0.00000 

0 Ck>000 

ESTIM. 

SIGMA (HCL) 

0.00839 

O OSG7 1 

0.00538 

0.00001 

0.00001 

O 00002 


parameter labels - X Y 2 XOOT YDOT 200T GKAV GRAV GRAV 

STATION/SATELLITE - O O O O O O O O O 

FILTER STATE (EPOCH) « 0.00330635 O.0O27rT-^5 O fK)28BSOI -0.0000024 1 0.00000423-0 00000009 0.00000000 -O.OOOQOOOQ 0.00000000 
FILTER SIGMAS (EPOCH)- 0.01444120 0.016526*10 p.O?!S3933 0.00002060 0.000017S1 0.00001337 0.00000000 0.00000000 0.00000000 


MBCAS I 

4 

0.00125162 

0.00401336 


VO 

45b 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


> MBIAS 2 REFRACT STA POS STA POS STA POS 

- 4 4 4 4 4 

• 0.00000011 0.00006394 '4^0001 1810 0.00015171 0.00018990 

- 0.00000099 0.00009t^T/>0. 00498390 0.00496892 0.00494967 


$$S$ RESIDUAL STATISTICS FOR TIME 
STATION 4 TYPE tO N- 

STATION 4 TYPE 13 N* 


3540. OOlf 
12 ME^N- 

12. MEt'^* 


3660.00 

.44700-06 SIGMA- 
.27950-03 SIGMA- 


. 12700-05 WEIGHTED SICMA- 
.11250-02 WEIGHTED SIGMA- 


.8684 

1.046 


FILTER OUTPUT AT TIME 3600.00 


NOMINAL TRAJECTORY (ECI )» 

4717.8187 

-2904.3301 

4392.7114 

-4.62769 

1 23220 

5.78336 

ESTIM. 

TOTAL (ECI) 

4717.8106 

-2904.3358 

4392 7041 

-4.62770 

1.23320 

5.76336 

EStIM. 

CORRECTION (HCL) « 

-0.00329 

-O 00302 

-0.01140 

-0.00000 

0.00001 

-O OOOOO 

ESTIM. 

SIGMA (HCL) 

0.00496 

0.00877 

0.00377 

o.ooorri 

0.00001 

0.00002 


parameter labels 

STATION/SATELLITE 
FILTER STATE (EPOCH) » 
FILTER SIGMAS (EPOCH)* 


X Y 2 YOOT YOOT 2001 GRAV GRAV 

O O O O O O O O 
0.00926432 0.00488644-0.00270946-0.00000097 O 00000750-0 0000035T 0.00000000-0 OOOOOOOO 
0.01290821 0.01590341 0.02307023 0.00002C83 O. 00001696 O 00001296 0.00000000 0.00000000 


GRAV MBIAS f 

O 4 

Q. OOOOOOOO O 005667S9 
0.00000000 0.00136934 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


• MBIAS 2 REFRACT ST.A POS STA POS STA POS 

• 4 4 4 4 4 

•-0 0000001I 0.000026490. 00029671 0-00037350 0.00036S44 

• 0.00000093 0.00007315 0.00497750 0.00495770 0.00494365 


O 

73 

Q 

Z 

» -0 


tr.%% RESIDUAL STATISTICS FOR TIME 
STATION 4 TYPE 10 N- 

STATION 4 TYPE 13 N- 


3660.00 TO 
12. MEAN* 
12. MEAN* 


3780.00 

. 13960-06 SIGMA* 
.39650-03 SIGMA- 


.10810-05 WEIGHTED SIGMA- .8436 

11920-02 WEIGIflEO SIGMA* 1.113 



m/SlIVESS AND TeaiN0L0GIC4L SYSTEMS, iNC 


FI'TIER OUlPUT AT TIME 3720.00 

trajectory (ECI)- 4125.7936 -2733.3490 5049.1735 
ESTiKi, TOTAL (ECI) - 4125 7945 -2733.3557 5049.1667 
£57:M. CORRECIION (HCU « >0.00334 -0.00397 >0.01213 
ESTIM. SIGMA (HCL) « 0.00425 O 00675 0.003B4 


-5.22598 

-5.22596 

- 0.00000 

3.00001 


1 6 1 358 5 14293 
I 61359 5.14293 
O 00001 -0.00001 
O OOOOl o oooot 


PARAMETER LABELS • X Y 2 XDOf 

STATION/SATELLITE » O O O O 

FILTER STATE (EPOCH) « O.OI0623B9 0.00554910-0.00243990-0 00000032 
FILTER SIGMAS (EPOCH)- 0.00847158 0.01395388 0.02186394 0.00002117 


VDOT ZOOT GRAV GRAV 

O O O O 
0.00000829-0.00000442 0.00000000 >0.00000000 
0.U0001695 0.00000756 0.00000000 0.00000000 


GRAV WIAS 1 

O 4 

0.00000000 0.00621968 
0.00000000 0.00080044 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FlLfER SIGMAS 


« MB I AS 2 REFRACT STA PQS STA POS STA POS 
-^4 4 4 4 4 

—0.00000071 O 00003685-0.00034657 O OOC42B68 0.00030179 
- 0.00000078 0.00006413 0.00497694 0.00495532 0.00491235 


%%%% RESIDUAL STATISTICS FOR TIME 
STATION 4 TYPE 10 N» 

STAllON 4 TYPE 13 N- 


3780.00 TO 3900.00 
II MEAN- .26350-06 SIGMA- 
It. MEAN- -.39110-03 SIGMA- 


.87 150-06 WEIGHTED SIGMA- 
.10040-02 WEIGHTED SIGMA- 


.8014 

.9356 


FILTER OUTPUT AT TIME 3840.00 


NOMINAL TRAJECTORY (EC1)=> 

3466 9873 

-2518. 1269 

5623.6961 

-5.73923 

1 96856 

4.41958 

ESTIM. 

TOTAL (ECI) 

3466.9770 

-2518. 1306 

5623.6885 

-5.73924 

1 96857 

4 41957 

ESTIM. 

CORRECTION (HCL) - 

-0.00080 

-0.00167 

-0.01326 

0.00000 

O OOOOl 

-0.00001 

ESTIM. 

SIGMA (HCL) 

0.00454 

0.00665 

0 00310 

0.00001 

O OOOOl 

0.00301 


iO 

tn 


parameter LABELS 
STATION/SATELLITE 
FILTER STATE (EPOCH) - 
FILTER SIGMAS (EPOCH)- 


X Y Z XDOT 

0 0 0 0 
0.01119814 0.00398565-0. 00366589 -0.0000008S 
0.00707678 0.01360886 0.02231765 0.00002163 


YOOT ZOOT CRAV GRAV GRAV MBIA5 I 

0 0 0 0 0 4 

0.00000747-0 00000547-0.00000000-0.00000000 0.00000000 0.00571648 
0.00001703 0.00000581 0.00000000 0.00000000 0.00000000 0.00062654 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


- I^IAS 2 REFRACT STA POS STA POS STA POS 

- 4 4 4 4 4 

—0.00000031 0.00003789-0.00030894 O.OOfT36995 0.00020514 
» 0.0000006? 0.0000S5B6 0.00497669 0.00495474 0.00490711 


%%%% residual STATISTICS FOR TIME 
STATION 4 TYPE 50 N* 

STATION 4 TYPE 13 N= 


3900.00 TO 4020.00 

3. MEAN- -.12890-05 SIGMA- .21000-05 WEIGHTED SIGMA- 1.455 

3. MEAN- .29920-03 SIGMA- .823f£^ <33 WEIGHTED SIGMA- 6638 


FILTER OUTPUT AT TIME 3960.00 


NOMINAL TRAJECTORY 

(ECI)* 

2752.1069 

-2362 1794 

6107 0271 

-6 15931 

2 29143 

3 62514 







ESTIM TOTAL (ECI) 

» 

2752.0951 

-2262.1790 

6107.0187 

-6. 15932 

2 29144 

3 62514 







ESTIM. CORRECTION 

(HCL) - 

0.00244 

Q 00151 

-0.01410 

O.OCOOO 

0 OOOOl 

-0 OOOOl 







ESTIM, SIGMA (HCL) 

• 

0.00488 

0 . 00728 

O 00438 

0.00001 

O OOOOl 

O OOOOl 







PARAMETER LABELS 

- X 

Y 

Z 

XDOT VDOT 

ZOOT GRAV 


GRAV 


GRAV 


MBIAS 

STATION/SATELLITE 

* 

0 

0 

O 

0 

0 

0 

0 


0 


0 



FILTER STATE (EPOCH) « 0.00747069 0.00299259-0.0038 1G94- 0.00000222 O 00000586-0.00000280 0.00000000 0.00000000 0.00000000 0.00515358 
FILTER SIGMAS (EPOCH)- 0.00664277 0.01362552 0.02291802 0.00002210 0.00001714 0.00000548 0.00000000 0.00000000 0.00000000 0.00054260 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


* MB I AS 2 REFRACT STA POS STA POS STA PQS 


- 4 4 4 4 4 

-- 0.00000024-0 00003459-0.00018796 O 00022862 0.00020SB8 

- 0.00000064 0.00003997 0.00497624 0.00495414 0.00490696 


ORIGINAL PAGE IS 
OF. POOR QUALITY 




AT 

TIME 

4080.00. 

PARAMETER 

lO. 

LABEL 

100004 

DELETED 

FROM 

STATE 

AT 

TIME 

4080 00. 

PARAMETER 

10. 

LABUL 

200004 

DELETED 

FROM 

STATE 

AT 

TIME 

4080.00. 

PARAMETER 

10, 

LABEL 

300004 

DELETED 

FROM 

STATE 

AT 

TIME 

4080,00. 

PARAMETER 

10. 

LABEL 

400004 

DELETED 

FROM 

STATE 

AT 

TIME 

4080.00. 

PARAMETER 

10. 

LABEL 

400004 

DELETED 

FROM 

STATE 

AT 

TIME 

4080.00. 

PARAMETER 

to. 

LABEL 

400004 

DELErtD 

FROM 

STATE 


FILTER CNJTPUr AT TIME 4080.00 

HOMINAL TRAJECTORY (ECI)« 1992.7473 -1999.6734 6491.4027 -6.47S58 
ESTIM. total CECI) « 1992.7342 -1969.6718 6491.3937 -6.47959 
ESTIM. CORRECTION (HCL) - 0.00480 0.00193 -0.01516 0.00001 
ESTIM. SIGMA (HCL.I « 0. 005 II 0.00829 0.00592 0.00001 


2.57706 2.77256 

2 . 57707 2 . 77255 
0-00001 - 0.00001 
O 00001 0.00001 


PARAMETER LABELS « X Y 2 XOOT YDOT 200r GRAV GRAV GRAV 

STATION/SATELLITE - OOOOOO OOO 

FILTER STATE (EPOCH) « 0.00746728 0.00299403-0.00380972-0.00000233 0.00000586-0.00000280 0.00000000 0.00000000 0.00000000 
FILTER SIGMAS (EPOCH)- 0.00664855 0.01365153 0.02351439 0.rXKK>3257 O 00001727 0.00000548 0.00000000 0.00000000 0.00000000 


FILTER OUTPUT AT TIME 4200.00 

NOMINAL TRAJECTORY (EC1)« 1201.2029 -1645.3574 6770.6657 -6.6S503 2 82090 1.87562 

ESTIM. TOTAL CECI) « 1201.1884 -1645.3545 6770.6555 -6.69504 2.82091 1.87561 

ESTIM. CORRECTION (HCL) « 0.00789 0.00213 -0.01598 0.00001 O 00001 -O.OOOOl 

ESTIM. SIGMA (HCL) - 0.00525 0.00949 0.00734 0.00001 O 00001 O OOOOl 

PARAMETER LABELS - X Y 2 XOOT VOOT 2t>OT GRAV GRAV GRAV 

STATION/SATELLITE - OOOOOOOO O 

filter STATE (EPOCH) - 0.00746450 0.00299476-0.00380600-0.00000223 O.OOOOOS86- 0.00000280 0.00000000 0.00000000 0.00000000 
FILTER SIGMAS (EPOCH)- 0.00665892 0.01367380 0.02409559 0.00002303 0,00001739 0.00000549 0.00000000 0.00000000 0.00000000 


FILTER OUTPUT AT TIME 4320.00 

NOMINAL TRAJECTORY (ECI)« 390 2692 -1294.4844 6940. 3561^ -6.80232 3.01909 0.94881 

ESTIM. TOTAL (ECI) - 390.2535 -1294.4801 6940.3448 -6.80233 3.01911 0.94879 

ESTIM. CORRECTION (HCL) - 0.01166 0.00203 -0.01625 0.00002 0.00001 -0.00001 

ESTIM. SIGMA (HCL) - O 00544 0.01081 0.00849 0.00001 0.00001 0.00001 

PARAME%'ER LABELS - X Y 2 XOOT YDOT 2D0T GRAV GRAV GRAV 

STATION/SATELLITE - OOOOOOOOO 

FILTER STATE (EPOCH) * 0.00746213 0.00299513-0.00380134-0.00000223 0.00000566-0.00000279 0.00000000 0.00000000 0.00000000 
FILTER SIGMAS (EPOCH)- 0.00667557 0.01369326 O 02465691 O 00002349 0.00001751 0.00000550 0.00000000 0.00000000 0.00000000 


FILTER OUrPUr AT TIME 4440.00 
NOMINAL TRAJECTORY (ECI)- -426.9611 

ESTIM. TOTAL (ECI) - -426.9777 

ESTIM CORRECTION (HCL) - 001S9B 

ESTIM. SIGMA (HCL) - 0.00583 


-922.7269 6997.7791 -6.79984 

-922.7210 6997,7653 -6.79985 

0.00160 -0.01567 0.00002 

0.01220 0.00926 0.00001 


3.16852 000699 

3.16853 000697 
0.00000 -O OOOOl 
0.00001 0.00001 


PARAMETER LABELS - X Y Z XOOT YDOT 200T GRAV GRAV GRAV 

STATION/SATELLITE -OOOOOOO OO 

FILTER STATE (EPOCH) » 0.00746004 0.0029 953 5 -0.0037 9656 -0.000002 24 0.000005B7-0.OD000279 0.00000000 0.00000000 0.00000000 
FILTER SIGMAS (EPOCH)- 0.00669897 O 01371054 0.02519511 0.00002394 0.00001764 O.OOOOOSS2 0.00000000 0.00000000 0.00000000 


[r 


BVSISESS AND TgCHN0U)GIC4L SYSTEMS, INC 


FILTER STATE tEPOCH) «-O.O1 196459 O.OOI5S85I 0. 01372915-0. CXXW1392 0.00000734 0.00001233-0. « X)00000 -0.00000000 0.00000000-0.00572779 
FILTER SIGMAS <EPOCH)» 0.00709953 0.02355629 0.16470665 0.00015902 0.00007499 0.00000475 OMjOOOOOO 0.00000000 0.00000000 0.00051519 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


« MBIAS 2 REFRACT STA POS STA POS STA P05 
31 31 31 31 31 

* 0.00000134-0.00009739 0.0G33707 2-0. 00 143922-0.00383002 
> 0.00000068 0.00004192 0.00389016 000481491 0.00419323 


%%%% RESIDUAL STATISTICS FOR TIME 
STATION 31 TYPE 10 N« 

STATION 31 TYPE 13 N- 


39660.00 TO 39780.00 

12- HEAN> 13900-06 SICHA« 

12. MEAN* .25090-03 SIGMA* 


.82170-06 WEIGHTED SIGMA* .6941 

.12660-02 WEIGHTED SIGMA* 1.199 


FILTER OUTPUT AT TIME 39720.00 


NOMiMAL TRAJECTORY 

(ECI)- 

1289.7009 - 

1668.0161 

6748.8733 

-6.69912 

2.75112 

1.96254 



E5TIM. TOTAL (ECI) 

* 

1289.7066 - 

1668.0199 

6748.8689 

-6.69912 

2 75111 

1 . 96256 



ESTIM. CGRRECTiaN 

(HCL) » 

0.00341 

-0.00637 

0.00379 

-0.00002 

o.oooco 

o ooooo 



ESTIH. SIGMA (HCL) 


0.00225 

0.00259 

0.00216 

0.00000 

0.00000 

o ooooo 



PARAMETER LABELS 

X 

Y 

Z 

XOOT 

(TOOT 2001 

GRAV 

GRAV GRAV 

MBIAS 1 

station/satellite 

“ 

O 

0 

O' 

0 

0 

0 

0 0 0 

31 


FILTER STATE (EPOCH) —0.01191212 0.00195711 O.OI 1 12740-0.00001062 0.00000559 O. 00001739 0.00000000-0.00000000 0.00000000-0.00539568 
FILTER SIGMAS (EPOCH)* 0.00724170 0.02345464 0.16543332 0.00015942 0.00007507 0-00000398 0.00000000 0.00000000 0.00000000 0.00036306 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


MBIAS 2 REFRACT STA POS STA POS STA POS 
31 31 31 31 31 

0.00000090-0. O0(X>9307 0 . 00256003-0.00 1 1 2262-0.00340266 
0.00000057 0.00(X)4095 0.00368938 0.00479503 0.00412870 


END OF MEASUREMENT DATA TAPE 


%%%% residual STATISTICS FOR TIME 
STATION 31 TYPE 10 N» 

STATION 31 TYPE 13 N= 


39780.00 TO 
7. MEAN- 
7 . MEAN* 


39900-00 

.56920-06 SIGMA- 
.46660-03 SIGMA* 


.11010-05 WEIGHTED SIGMA* 
.15010-02 WEIGHTED SIGMA- 


.9007 
1 . 365 


FILTER OUTPUT AT TIME 39840-00 


NOMINAL TRAJECTORY 

(EC1)» 

477.5626 

-1325.3116 

6929.1433 

-6.81828 

2.95292 

1.03798 

ESTIM. 

TOTAL (ECI) 

« 

477.5681 

-1325.3157 

6929. 1406 

-6.81828 

2 95291 

1 .03799 

ESTIM- 

CORRECTION 

(HCL) - 

0.00114 

-0.00575 

0.00442 

-0.00002 

0 CKKXX) 

0.00000 

ESTIM. 

SIGMA (HCL) 

» 

0.00220 

0.00265 

0.00228 

0.00000 

0 OOOOO 

O OOOOO 


o 2 

m 50 
•tj Q 


o 2 
o > 

50 r- 

o -o 

c; > 
> Q 
r- m 


-< c /> 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE (EPOCH) 
FILTER SIGMAS (EPOCH) 


* X Y Z XOOT YOOT ZOOT GRAV GRAV 

*00000000 
*-0.01028879-0.00360944 0.04112749-0.00004005 0.00001975 0.00000865-0.00000000-0.00000000 
= 0.00739320 0.02360592 0.16738533 0.00016126 0.00007591 0.00000342 0.00000000 0.00000000 


GRAV MBIAS I 

O 31 

0.00000000-0.00524538 
0.00000000 0.00032576 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
FILTER SIGMAS 


MBIAS 2 REFRACT STA POS STA POS STA POS 
31 31 31 31 31 

0.00000052-0.00005719 0.00221876-0.00090583-0.00338876 
0.00000052 0.00002683 0.00363468 0.00478666 0.00412166 


SUMMARY OF RESIDUAL STATISTICS FOR 1082. MEASUREMENTS 

WEIGHTED SIGMA * 1.0160 . SUM(RES*«2/RVAR) * 1116.8 , SUM( ALOGiRVAR )) * -21929. 


ALOG LIKELIHOOD « 10406- 
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i 


SMOOrHER OUTPUT AT TIME 39730.00 


PARAMETER LABELS « X V 2 XOOT 

STATION/SATELLITE * O O O O 

FILTER STATE fEPOCH) >'*0.01191213 0.00195711 0.01 1 1 3740^0. 00001062 
SMOOTHER ST ATE (EP0CH)«~O. 01039023-0. 00384667 0.04291335-0.00004176 
SMOOTHER SIGMAS(EPOCH) 0.00719136 0.03333254 0.16449102 0.00015650 


VOOT 200T GRAV GRAV GRAV MBT^S 1 

O O O O O 31 

0.00000559 0.00001239 0.00000000-0.00000000 0.00000000-0.00539568 
0.00002055 0.00000865-0.00000000-0.00000000 0.00000000-0.00524538 
0.00007463 0.00000339 0.00000000 0.00000000 0.00000000 0.00032576 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


MBIAS 2 REFRACT STA POS STA POS STA POS 
31 31 31 31 31 

0.00000090-0.00009307 0.00256003-0.001 12263-0.00340266 
0.00000052-0.00005719 0.0022 IB76-0.00090563-0. 00338876 
0.00000052 0.00002683 0.00363468 0.00478666 0.00412166 


NOMINAL TRAJECTORY (ECI|- 
SMOOTHER TOTAL (ECll 
FILTER ESTIM. fHCL) 
SMOOTHER ESTIM. CHCL) 
SMOOTHER SIGMA (HCL) 


1389.7009 
S2B9.7065 
0.0034 109 
0.0036475 
0.0022143 


-1668.0161 

-1668.0197 

-0.0063667 

-0.0061407 

0.0025913 


6748.8733 

6748.8687 

0.0037916 

0.0036228 

0.0021084 


-6.69912 

-6.69912 

-0.0000196 

-0.0000162 

0.0000040 


2.751 12 
2.75III 
O.OOOOQT7 
0.0000013 
0.0000030 


1 . 96254 
1.962S6 
0.0000045 
0.GO00032 
0.CO00025 


SMOOTHER OUTPUT AT TIME 39200.00 


to 

00 


PARAMETER LABELS - X Y 2 XDOT YDQT 200T GRAV GRAV GRAV MBIAS 1 

STATION/SATELLITE - O O O O O O O O 0 31 

FILTER state (epoch) —0.01196459 0.00155251 G. 01372915-0. 00001392 0.00000734 0.00001233-0.00000000-0.00000000 0.00000000-0.00572779 
SMOOTHER STATE(EPaCH)«-0. 01048145-0. 00404558 0.04442501-0.00004330 0.00003123 0 . 00000866 - 0 . 00000000 - 0.00000000 0 . 00000000 - 0.00524538 

SMOOTHER SIGMAS(EPOCH) 0.00699729 0.03'.f83163 0.16142531 0.00015558 0.00007.127 0.00000337 0.00000000 0.00000000 O.QOOOOOOO 0.00032576 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


MBIAS 2 REFRACT STA POS STA POS ST.^i POS 
31 31 31 31 31 

0.00000134-0.00009739 0.00337072-0.00143932-0,00383002 
0.00000053-0.000057 19 0.00321876-0.00090583-0.00338876 
0.00000052 0.00002683 0.00363468 0.00478666 0.00412166 


NOMINAL TRAJECTORY (EC1)> 
SMOOTHER TOTAL (ECI ) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


2081.0206 

2081.0262 

0.0068344 

0.0060435 

0.0034619 


-1983.7978 

-1983.8006 

-0.0065368 

-0.0062858 

O.OC25332 


6459.3856 

6459.3791 

0.0014804 

0.0021863 

0.0018484 


-6.47175 

-6.47175 

-0.0000193 

-0.0000161 

0.0000C(39 


2.50479 
2 . 50479 
0.0000007 
- 0.0000000 
0-0000030 


2.85583 

3.85584 
O 0000069 
O 0000061 
O 0000038 


SMOOTHER OUTPUT AT TIME 39480.00 

PARAMETER LABELS « X Y 2 XOOT YDOT ZOOl GRAV GRAV GRAV MBIAS I 

STATION/SATELLITE > O O O O O O O O O 31 

FILTER STATE (EPOCH) >-0.01230410 0.00525684-0 . OOB26364 0.00000937-0.00000419 O.OOOOI3J9 0.00000000 0.00000000-0,00000 0 00-0,00506149 
SMOOTHER ST ATE( EPOCH )»-0. 01056495-0. 00422756 0.04581033-0.00004452 0.00002183 0.00000867-0.00000000-0.00000000 0.00000000-0.00524538 
SMOOTHER SIGMASlEPQCH) 0.00680871 0.02241107 0.15823161 0.00015255 0.00007185 0.00000336 0.00000000 0.00000000 0.00000000 0.00032576 

PARAMETER LABELS 
STATlOfZ^SATELLlTE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 

NOMINAL TRAJECTORY (ECI )= 2838 7395 -2267,5407 6065.3039 -6.13967 2.21784 3.70344 

SMOOTHER TOTAL (ECI) « 3838.7348 -2267.5428 6065.2944 -6.13966 2.21784 3.70345 

FILTER ESTIM. (HCL) > 0.0081388 -0.0064660 0.0006029 -0.0000192 -0.0000025 0.0000130 


» MBIAS 2 REFRACT STA POS STA POS STA POS 
= 31 31 31 91 31 

> 0.00000017-0.00010301 0. 00162578-0. 00079WS 1-0.002867 10 

> 0.00000052-0.00005719 0.00221876-0.00090583-0.00338876 
0.00000052 0.00002683 0.00363468 0.00478666 O.0O412166 


o o 

TI 2 

“0 Q 
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33 r- 
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SMOOTHER ESTIM. CHCL) 
SMOOIMR sigma IHCLI 


» O 0091972 -O G06I1S6 
- O 0039091 O 0024539 


0.0001315 O 0000155 >0 0000017 O 000009 1 
O 0015979 O 0000037 O 0000029 O 0000032 


SMOOTHER OUTPUT AT TIME 39390 00 

PARAMETER LABELS -X Y 2 XOOT 

STATION/SATELLITE - O O O O 

FILTER STATE (EPOCH) >'0.014 10C99 000027999 O 03ai34G5'0 000037 10 
SMOOTHER ST ATE ( EPOCH) -'0.0 106354 2'0. 004 37357 004693099'0 00004559 
SMOOTHER SIGMAS(EPOCH) 0.00692439 O 03197792 0.15495235 0.00014944 


yOOT ZOOT GRAV GRAV GRAV MIAS f 

O O G O O 31 

0.00001290 0.00001394 O.OOOOOOOO'O OOOOOOOO'O OOOOOOOO'0 . 00017022 
0.00002233 0. 00000969-0. OOOOOOOO'O OOOOOOOO O OOOOOOOO'O 00534539 
O 00007040 0.00000339 O 00000000 O 00000000 O OCOOOOOO O 00032576 


PARAMETER LABELS 
STATION/SATELLiTE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


MB I AS 2 REFRACT St A PQS ST A POS STA POS 

31 31 31 31 31 

O . OOOOOOOO'O 00019162 0.001 35949 -O 00094 712-0 00243527 
0. 00000052-0. 000057 19 O 00221879-0 00090583-0 00338876 
0.00000052 O 00002693 O 00363469 O 00479966 O <^12169 


NOMINAL TRAJECTORY (ECD- 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM (HCL) 
SMOOTHER ESTIM (HCL) 
SMOOTHER SIGMA (HCL) 


35505579-25149399 55779393 -570907 169493 
35505939-2514.6390 55729279 -570907 199492 
0.0091577 '0.0025187 -O 0049125 -0.0000160 -O 0000025 
0.0099593 '0.0055595 -0.0024915 -O 0000141 -C 0000039 
0.0031416 0.0023393 0 0015530 0.0000034 0 0000029 


4 49164 
4 49165 
OOG 00130 
O 00)0121 
O 000036 


SMOOTHER OUTPUT AT TIME 39240 OO 


parameter labels « X Y Z XCOT YDCT ZOO I GRAV GRAV GRAV 

STATIOM/SATELLITE - O O O O O O O O O 

^ FILTER state (EPOCH) —0.01265392 0.00701390-0 01926674 0.0000 1396 -0.0000067 2 0) 00001497-0 00000000-0 OOOOOOOO'O OOOOOOOO 
S SMOOTHER 3TATE(EP0CH)--0 01067324-0. 00442226 O 04734736-0 00004596 0.00002250 0* 00000970-0 00000000-0 00000000-0 OOOOOOOO 
SMOOfHER SXGMAS(EPOCH) 0.00644052 0.02155520 O 15172493 0.00014639 0.00009999 C‘ O0C0O337 0.00000000 O OOOOOOOO O OOOOOOOO 


NOMINAL TRAJECTORY (EC1)» 
SMOOTHER TOTAL (ECt) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM. (HCLl • 

SMOOTHER sigma (HCL) « 


4304 9572 -2721 0619 4990 2020 -5.1B379 1 54090 
4204 9617 -2721 0621 4990 189S -5.18379 1 54099 
0.0292320 'O 0139232 -O 010I673 -0.0000231 -O 0000122 
O 0111706 -O 0045501 -O 0055969 -0.0000120 -0 0000063 
O 0033816 O 0021911 0 00)9714 0.0000030 O 00CGO2O 


5 20759 
5 30759 
O 0000339 
O 0000147 
O 0000040 


SMOOTHER OUTPUT AT TIME 39130 OO 


OO 

T| 2 

O 

o 

lO ^ 

d !>' 

£I ^ 
< u. 


PARAMETER LABELS > X Y Z XOOT YOOT 2D0T CRAV GRAV GRAV 

STATION/SATELLITE » O O O O O O O O O 

FILTER STATE (EPOCH) *-0 01265457 O 0070 1294-0 01625917 O 00003395-0 00000672 O 00001497-0 00000000*0 OOOOOOOO'O OOOOOOOO 

SMOOTHER STATC(EP0CM)—0 01069227-0 00439753 O 047 1595 1-0 00004577 O 00002241 0 0CKXX)ft72 O 00000000-0 OC)000000-0 OOOOOOOO 

SMOOTHER SIGMAS(EPOCH) O 00625964 0.02114227 O 14955694 O 00014341 O 00006759 O 00000340 O OOOOOOOO O OOOOOOOO O OOOOOOOO 


TOMINAL TRAJECTORY (ECI)> 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM. (HCL) 
SAOOTHER ESTIM (HCL) 
SMOOTHER SIGMA (HCL) 


4791 2909 
4791.2945 
O 0291601 
O 0116737 
O 0034797 


-2993 4429 
-2893 4422 
-O 0107151 
'O 0030324 
O 0019970 


4326 . 4909 
4326 4663 
-O 0194131 
-00090039 
O 0025013 


-4 57516 
-4 57515 
-O 0000157 
-O 0000090 
O 0000026 


I 16 ITS 
1 161T4 
'O 0000194 
-0 0000092 
O 0000032 


S 93953 
5 83955 
O 0000373 
O 0000168 
O 0000044 


o 


SMOOTHER OUTPUT AT TIME 

39000 00 









PARXMETER LABELS 

X 

Y 

Z 

XOOT 

YOOT 

ZOOT 

CRAV 

GBtV 

CRAV 


vsi 'fimiisis 'iroiofrionimj, asnmmw 


StIOOTHER OUTPSJT AT TIME 34320.00 



PARAMETER LABELS * X V Z XOOT VDOT ZDOT GRAV GRAV GRAV 

STATION/SATELLITE *^000000000 
PILTER STATE (EPOCH) —0.01288091 0.00667978-0.01551303 0.00001 107-0 00000539 O 00001506-0.00000000-0.00000000-0.00000000 
SMOOTHER STATE(EP0CH)--O. 01392362-0. 00775946 0.07764766-0.00007636 O 00003731 00(X}OI060'0. 00000000 0.00000000-0.00000000 
SMOOTHER SIGHAS(EPOCH) 0.00493206 O 01242826 O 08935511 0.00008583 0.00004028 0.00000221 O OOOOOOOO 0.00000000 0.00000000 


NOMINAL TRAJECTORY (ECD- 
SMOOTHER TOTAL (ECI ) 
FILTER ESTIM. (HCL) 
SMGOTHEn ESTIM. (HCL) 
SMOOTHER SIGMA (HCL> 


-2223.5573 -52.8471 

-2223.5492 -52.8539 

-0.0089958 -0.0045663 
-0.0068559 -0.0053091 
0.0031333 0.0028263 


6713.9066 -6.41486 
6713.9084 -6.4 1486 
0.0057511 -0.0000201 
0.0062154 -0.0000159 
0.0014741 0.0000037 


3 27866 -2 09260 

3.27666 -2.09258 

O 0000098 -0.00'0152 
0.0000082 -o.oo::oi34 
0.0000026 0.0000015 


SMOOTHER OUTPUT AT TIME 34200.00 

PARAMETER LABELS > X Y Z XOOT YOOT ZOOT GRAV GRAV GRAV 

STATION/SATELLITE « O O O O O O O O O 

FILTER STATE (EPOCH) —0.01292766 0.00666S88-0. 01531427 0.00001087-0.00000530 0.000015 10-O.OOOOOCK)0-0.00000(X>0-O.OOOOOCOO 
SMOOTHER STATE(EPOCH)— 0.01410321-0.007626.^7 0.07735176-0 0(XX)78 15 0.00003717 0.00001082-0.00000000 0.00000000-0.00000000 
SMOOTHER SIGMAS(EPOCH) 0.00480205 0.01244604 0.08008353 0.(XXX>B555 0.00004017 O.00G0O222 0.00000000 0.00000000 0.00000000 


NOMINAL TRAJECTORY (ECI 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM, (HCL) 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


- 1437 . 9 158 -444 . 5628 

-1437,5068 -444.5693 

-0.0055320 -0.0062105 
-0.(X>38a83 -0.0067623 
0.(X)27318 0.0026839 


6910.0647 -6.66154 
6910.0639 -6.66155 
0.0078919 -0.0000227 
0.0079470 -0.0000188 
0.0014850 0.0000034 


3.24315 -1.17230 

3.24314 -1.17227 

0.0000093 -0,0000112 
0.0000081 -0.0000104 
0.0000027 0.0000078 


g SMOOTHER OUTPUT AT TIME 34080.00 


o o 

-TI g 

*0 52 

O 2 
O > 
r! 


PARAMETER LABELS - X Y Z XDOT YDOT ZOOT GRAV GRAV GRAV 

STATION/SATELLITE - OOOOOO OO O 

FILTE;: state (epoch) »-0.012S7605 0.00667141-0.01518205 0.00001073-0.00000524 O. OOOO IS 15-0. OOOOOOOO -0.00000000-0. (X)000000 
SMOOTHER STATE! EPOCH)«-0.01<325524-0 00742048 0.07648873-0.00007737 0.00003677 O 00001105-0.00000000 0.00000000-0.00000000 
SMOOTHER SIGMAS! EPOCH) 0.00467897 0.01246169 0.08881470 0.00008530 0.00004007 0.00000223 0.00000000 0.00000000 0.00000000 


NOMINAL TRAJECTORY (ECI )- 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM, (HCL) 
SMOaTHER SIGMA (HCL) 


'629.0865 

-629.0767 

-0.0015608 

-0.0003099 

0.0023868 


-829. 1093 
-829.1155 
-0.0078040 
-0.0081892 
0.0025538 


6994.4987 
6994.4953 
O 0090308 
0.0088358 
0,0014495 


-6.80060 

-6.80081 

-0.0000250 

-0.0000215 

0.0000031 


3 15733 
3, 15733 
0.0000086 
0.0000076 
0.0000027 


-O 23302 
-0.23300 
-O 0000070 
-0.0000069 
0.0000022 


«o -n 

g: > 

C ^ 

j — 
-< <flr 


i 

I 

■i 


PARAMETER RECONSTRUCTION, NO ^ 9 N ^ 15 

OLD LABELS > 123456 

NEW LABELS > 123456 


8 8 8 

888 100031 200031 300031 4PCC3I 400031 400031 


SMOOTHER OUTPUT AT TIME 33960.00 


PARAMETER LABELS ^ X Y Z XOOT YOOT 

STATION/SATELLITE « O O O O O 

FILTER STATE (EPOCH) *-0.01302252 0,00669649-0.01520233 O.COOO 1074-0.00000526 
SMOOTHER STATE! EP0CH)=-0.01436a 12-0.007 10705 0.07485697-0.00007583 0.(XX>03€01 
SMOOTHER SIGMAS(EPOCH) O.O04S6BOI 0.01247485 0.08856004 0.0000BS07 0.00003997 


ZOOT GRAV GRAV GRAV MB I AS 1 

O O O O 31 

O . OOOO 1 520-0 . 00000000-0 . 00000000-0 . 00000000-0 . 005 16 I 55 
O 00001128-0 00000000-0 00000000-0.00000000-0.00505704 
0.00000226 0.00000000 0.00000000 0.00000000 0.00047105 


PAI3AMETER LABELS 
STATION/SATELLITE 
FILTER STATE 


MBIAS 2 REFRACT STA POS STA POS STA POS 
31 31 31 31 31 

0.00000097 >0.00013497-0.006S7906 0.00014390-0.00368113 
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SMOOTHER STATE 
SMOOTHER SIGMAS 


- 0,00000058-0.00012990-0.00740505 0.00008527-0.00344118 
0.00000050 0.00003847 0.CO4II089 0.00499719 0.00338917 


NOMINAL TRAdECTORV (ECI| 
SMOOTHER TOTAL (ECt) 
FILTER ESTIM. (HCLl 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


189.8874 

189.8978 

0.0027560 

0.0035906 

0.0021466 


-1200.2840 

-1200.2898 

-0.0092S31 

-0.0094970 

0.0024374 


6965.8335 

6965.8272 

0.0090843 

0.0087669 

0.0013534 


-6.83039 

-6.83039 

-0.0000268 

-0.0000238 

0.0000028 


3.02058 

3.02057 

0.0000076 

0.0000068 

0,0000028 


0.71017 
0.71020 
-O 0000024 
-O 0000027 
O 0000025 


SMOOTHER OUTPUT A'i TIME 33840.00 


PARAMETER LABELS « X V Z XOOT 

STATION/SATELLITE - O O O O 

FILTER STATE (EPOCH) «-0. 01025805-0. 00352851 O. 03989700-0. OOOC3869 
SMOOTHER STATE(EP0CH)«-0. 01443920-0.00670008 0.07251688-0.00007359 
SMOOTHER S)GMAS(EPOCH) 0.00447487 0.01248474 0.08833056 0.00009487 


Yoor zoor gpav grav grav mbias i 

O O O O O 31 

0.00002027 0-000008 79- 0.00000000 0.00000000 0.00000000-0.00420297 
0.00003493 0.00001 f 50-0.00000000 0.00CKXXXX)-0. 00000000-0.00505704 
0.00003989 O 00000230 0.00000000 0.00000000 0.00000000 0.00047105 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


- MBIAS 2 REFRACT STA POS STA POS STA POS 

- 31 31 31 31 31 

-0.00000099 0.00008011-0.00238554 O 00012084-0.001 I 1765 
« 0.00000058-0. 00012990-0. 0074053S 0.00008527-0.00344118 

0.00000050 0.00003847 0.00411089 0.00499719 0.00338917 


NOMINAL TRAJECTORY (ECI )« 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


1005. 7S64 
1005.8072 
0.0035039 
0.0076911 
G.C020379 


-1552.0361 

-1552.1013 

-0.0096903 

-0.0105755 

0.0023340 


6824,5052 

6824.4959 

0.0069416 

0.0076408 

0.0012256 


-6.74978 

-6.74978 

-0.0000208 

-0.0000255 

0.0000027 


2 83505 
2 83504 
0.0000023 
0.0000055 
O 0000029 


I 64218 
1.64221 
O 0000033 
O 0000020 
0.0000026 


SMOOTHER OUTPUT AT TIME 33720.00 


0 2 

•n 2 

O S 
O > 
73 n 


parameter labels -X Y z XOOT 

STATION/SATELLITE - 0 O O O 

FILTER STATE (EPOCH) --0.01 173550-0.00031796 0.02422561-0 00002062 
SMOOTHER STATE(EPOCH)— 0.01454938-0.00648348 0.07147197-0.00007257 
SMOOTHER SIGMAS(EPOCH) 0.00439642 0.01248467 0.08806834 0.00008464 


YDOT ZOOT GRAV GRAV GRAV MBIAS 1 

O O O O O 31 

0.00001225 0.00001148 0.00000000-0.00000000-0. 00000000-0. 00317798 
0.00003443 0.00001169-0.00000000-0.00000000-0.00000000-0.00505704 
0.00003980 0.00000234 0.00000000 0.00000000 0.00000000 0.00047105 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


« MBIAS 2 REFRACT STA POS STA POS STA POS 
- 31 31 31 31 31 

— O.0Ck)O0O33 0.00007278-0 00029536-0.00002773 0.00050737 
« 0. 00000058 -O.(X)012990-0. 00740535 0.00008527-0.00344119 
0.00000050 0.00003847 0.00411089 0.00499719 0.00338917 


lO :g 

<= 5 
a 

C ^ 

^ 03 


NOMINAL TRAJECTORY (ECI)» 
SMOOTHER TOTAL (ECI) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


1805.4718 -1878.8608 
1805 4829 -1878 8654 
O 0041293 -0.0097481 
0.0117773 -0.0112952 
0.0020560 O . 0022404 


wri 7555 
6572.7431 
0.0052776 
O . 0053870 
O 0011324 


-6.56017 

-6.56017 

-0,0000215 

-0.0000264 

0.0000028 


2 60367 
2.60367 
-O 0000013 
0.0000035 
0.0000030 


:2, 54806 
2.S4809 
O 0000094 
O 0000074 
O 0000026 


SMOOTHER OUTPUT AT TI»E 33600.00 


PARAMETER LABELS « X Y Z XOOT YOOT 

STATION/SATELLITE * 0 0 O O O 

FILTER STATE (EPOCH) »-0. 0105426 1-0. 00333562 0.04032661-0 00003641 0.00001992 
SMOOTHER SfATE(EPQCH)>-0 0I4646B9-0.0062S672 0.07034010-0.00007145 0.00003387 
SMOOTHER SIGMASC EPOCH) 0.00433628 0.01246784 0.08770679 0.00008433 0.(XX)03966 


2D0T GRAV GRAV GRAV MBIAS I 

O O O O 31 

0.00000948-0. 00000000-0. (XXXXXXX) 0.00000000-0.0032Ci49 
O . 00001 187 -O , 00000000-0.00000000-0 00000000-0 , 005057G4 
0.00000239 0.00000000 0.00000000 O 00000000 0.00047105 
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PARAMETER LABELS - MBIAS 2 REFRACT STA POS STA POS STA POS 
STATION/SATELLITE » 31 31 31 31 31 

FILTER STATE « 0.00000001 0.00010116-0.00022424 0.00001401 0.0004? t2 I 

SMOOTHER STATE • 0.0000005B- 0.000 1 2990 -0.00740523 O 00009527-0.00344 1 1B 

SMOOTHER SIGMAS 0.00000050 0.00003B47 0.00411009 0.00499719 0.0033B917 

NOMINAL TRAJECTORY f ECl )« 257S.9927 -2175.2809 6214.6000 -6.26450 2.330lfl 3.41322 

SMOOTHER TOTAL (ECI) = 2S76.O037 -2175.2926 6214.5843 --6.26449 233010 3.41325 

FILTER ESTIM. (HCL) * 0.0071108 -0.0094902 0.0020559 -0.0000104 -0.0000036 0.0000115 

SMOOTHER ESTIM. CHCL) » 0.0155920 -0.0115129 0.0019B52 -0.0000262 O.OOOOCOft 0.0000131 

SMOOTHER SIGMA (HCL) * 0.0021656 0.0021503 0.0011745 0.0000030 0.0000029 O 0000027 


SMOOTHER OUTPUT AT TIME 33400.00 


PARAMETER LABELS - X Y 2 XOOT YOOT 2D0T GRAV GRAV GRAV MBIAS 1 

STATION/SATELLITE « O G O O O O O O O 31 

FILTER STATE (EPOCH) *-0.01185935-0.00676767 0.06460779-0.00006069 0.00003149 O 00000923-0. C3000000-0. 00000006 0.00000000 0.00015798 
SMOOTHER ST ATE( EPOCH) *-0.0 1469780-0. 005932 18 0.06844481-0.00006959 0.00003298 O 00001204-0.00000000 0.00000000-0.00000000-0.00505704 
SMOOTHER SIGMAS(EPOCH) 0.00425662 0.01243498 0.08722282 0.00008390 0.00003947 0.00000245 0.00000000 0.00000000 0.00000000 0.00047105 


PARAMETER LABELS 
STATION/SATELLITE 
FILTER STATE 
SMOOTHER STATE 
SMOOTHER SIGMAS 


- MBIAS 2 REFRACT STA POS STA POS STA POS 
• ' 31 31 31 31 31 

—0.00000002-0.0000^51 0.00020439 0.00002857 0.00137975 
> 0.00000058-0.00012990-0 00740535 0.00008527-0.00344118 
0.00000050 0.00003547 0.00411009 0.00499719 0.00338917 


NOMINAL TRAJECTORY (ECI)« 3304.0902 -2436.5711 5755.7709 -5.06736 2 0!B6« 4 22368 
SMOOTHER TOTAL (ECI) * 3304.9008 -2436.5736 5755.7518 -5.86735 2.01867 4.22371 
FILTER ESTIM. (HCL) « 0.0084830 -0.0069912 -0.0018416 -0.0000154 -O 0000047 0.0000115 
SMOOTHER ESTIM. (HCL) « 0.0188414 -0.0110819 -0.002509? -0.0000246 -0.0000027 O 0000188 
SMOOTHER SIGMA (HCL) * 0.0023157 0.0020543 0.0014257 0.0000031 O 000002G O 00C0029 


O 

ro 


SMOOTHER OUTPUT AT TIME 33360.00 

PARAMETER LABELS * X Y 7 XOOT YOOl 20OT GRAV GRAV GRAV 

STATION/SATELLITE * 000000000 

FILTER STATE (EPOCH) —0.00809563 0.00226557-0.00269560 0.00000623 0.00000016 G 00000884 O.ODOOOOOO 0.00000000 0.00000000 
SMOOTHER STATE(EP0CH)«-0. 01478549-0.00580487 0.06788934-0.00006899 0.00003268 0.000012 17-0. OOOCOOOO 0.00000000-0.00000000 
SMOOTHER SIGMAS(EPOCH) 0.00418174 0.01238971 0.08662333 0.00008338 0.00003924 0.00000252 0.00000000 0.00000000 0.00000000 

NOMINAL TRAJECTORY (ECI)- 3980.3483 -2SS8. 4563 6203.6326 -5.37501 1 67436 4.96626 

SMOOTHER TOTAL (ECI) * 3980.3581 -2658.4575 5203.6099 -5.37501 1.67435 4.96629 

FILTER ESTIM. (HCL) » 0.0175172 -0.0141730 -0.0011385 -O 0000156 -O 0000136 0.0000245 

SMOOTHER ESTIM. (HCL) * 0.C2I2I09 -0.0098631 -0.0079489 -0.0000213 -O 0000071 O 0000240 

SMOOTHER SIGMA (HCL) * 0.0024536 0.0019452 0.0018710 0.0000031 0.000002T O 0C00032 


SMOOTHER OUTPUT AT TIME 33240.00 

PARAMETER LABELS » X Y Z XOOT YOOT 200T GRAV GRAV GRAV 

5TATI0N/SATELLEVE « O O O O O O O O O 

FILTER STATE (EPOCH) —0.00809529 0.00226591-0.00269871 0.00000623 0.00000016 0,00000884 0.00000000 0.00000000 0.00000000 
SMOOTHER STATE(EPOCH)— 0.01489570-0.00581037 0.06823271-0.00006924 0.00003278 0,00001227-0.00000000 0.00000000 -O.DOOOOOOO 
SMOOTHER SIGMAS(EPOCH) 0.00410578 0.01233420 0.08594364 0.(XXX>8279 0,0CO0389B 0.00000260 0.00000000 0.00000000 0.00000000 

NOMINAL TRAJECTORY (ECU* 4591.3958 -2637.3228 4567.0704 -4.79528 1.30267 5 62882 

SMOOTHER TOTAL (ECI) * 4591.4044 -2837.3225 4567.0442 -4.79526 1 30265 5.62885 

FILTER ESTIM. (HCL) » O 0IB8S32 -0.0117729 -0.0069333 -0.0000103 -O 0000173 O 0000255 
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SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA <HCL) 


» 0.0223922 -0.0077423 -0.0140761 -0.0000164 -0.0000121 0.000029^3 
- 0.002S32S 0.0016270 0.0024499 0.0000029 0.0000028 0.0000036 


SMOOTHER OUTPUT AT TIME 33120.00 


PARAMETER LABELS « X Y Z XDOT VDOT ZDOT GRAV GRAV GRAV 

STATION/SATELLITE >> OOOOOOOOO 

FILTER STATE (EPOCH) --0.0060949S 0.00226631-0.00270220 0.00000623 0.0000001S 0.00000664 0.00000000 0.00000000 0.00000000 
SMOOTHER STATED EP0CH)«-0. 01501208-0.00569059 0.06906261-0.00006995 0.00003310 O.O00C123S-0.O0O0OO0O 0.00000000-0.00000000 
SMOOTHER SIGMAS(EPOCH) 0.00403126 0.01226666 0.06520073 0.00006214 0.00003EC9 0.00000267 10.00000000 0.00000000 0.00000000 


NOMINAL TRAJECTORY (ECI)- 
SMOOTHER TOTAL (ECl) 
FILTER ESTIM. (HCL) 
SMOOTHER ESTIM. (HCL) 
SMOOTHER SIGMA (HCL) 


5136.0676-2970.2404 3856.3545 -4.13741 0.90958 6.20046 
5126.0944 -2970.2364 3856.3246 -4.13740 0.90957 6.20048 
0.0186020 -0.0086115 -0.0126916 -0.0000045 -0.00002C9 0.0000251 
0.0221167 -0.0046494 -0.0205240 -0.0000099 -0.0000175 0.0000307 
0.0025162 0.0017293 0.0030997 0.0000025 0.0000029 0.0000040 


SMOOTHER OUTPUT AT TIME 33000.00 

parameter LABELS - X Y Z XDOT YCOT ZOOT GRAV GRAV GRAV 

STATION/SATELLITE *0 0 0000000 

FILTER STATE (EPOCH) *-0.00609462 0.00226676-0.00270609 0.00000624 0.00000015 0.00000664 0.00000000 0.00000000 0.00000000 
SMOOTHER STAT£( EPOCH)— 0.01512492-0.00601130 0.07013711-0.00007067 0.00003353 0.00001241-0.00000000 0.00000000-0.00000000 
SMOOTHER SIGMAS(EPOCH) 0.00395969 0.01219362 0.06440394 0.00006146 0.00003836 0.00000273 0.00000000 0.00000000 0.00000000 

NOMINAL TRAJECTORY (ECI)- 5581.6721 -3055.0214 3082.9785 -3.41202 O 50148 6 67171 

SMOOTHER TOTAL (ECI) - 5581.6766 -3055.0174 3062.9454 -3.41200 0.50146 6.67174 

FILTER ESTIM. (HCL) " 0.0173428 -0.0047515 -0.0180432 0.0000014 -0.00002/7 0.0000232 

SMOOTHER ESTIM. (HCL) * 0.0201933 -0.0005754 -0.0268249 -0.0000023 -0.000033V O 0000309 

SMOOTHER SIGMA (HCL) * O 0023996 0.0017192 0.0037561 0.0000022 O 0000033 0.0000041 


SMOOTHER 0U1PUT AT TIME 32880.00 

PARAMETER LABELS « X Y Z XDOT YDOf ZOOT GRAV GRAV GRAV 

station/satellite * OOOOOOO O O 

FILTER STATE (EPOCH) «-0. 00809429 0.00226731-0.00271041 0.00000624 0.00000015 O.OOOtJS>1;£:4| 0.00000000 0.00000000 0.00000000 
SMOOTHER STATE(EP0CH)«-0. 01522786-0. 0061512B O 07130463-0.00007 166 0.00003400 0.00001245-0.00000000 0.00000000-0.00000000 
SMOOTHER SiGMAS(EPOCH) 0.00369236 0.01210838 0.06356013 0.00008074 0.00003805 0.00000279 C.OOOOOOOO 0.00000000 O.OOOOCXIOa 

NOMINAL TRAJECTORY (ECI)* 5944.7403 -3090.3596 2259.4758 -2. 63086 G 08500 7 03474 

SMOOTHER TOTAL (ECI) * 5944.7419 -3090 3535 2359.4397 -3.63083 O.OB498 7.03476 

FILTER ESTIM. (HCL) « O 0145296 -0.0002952 -0.0326102 0.0000072 -0.0000271 O 3000197 

SMOOTHER ESTIM. (HCL) * 0.0165448 0.00441.16 -0.03243B7 0,0000059 -O. 0000283 O 00002 85 

SMOOTHER SIGMA (HCL) * 0.0022083 0.0018824 0.0043501 0.0000022 0.0000038 0.0000040 


SMOOTHER OUTPUT AT TIME 32760.00 

PARAMETER LABELS - X Y Z XDOT YDOT ZDOT GRAV GRAV MV 

STATION/SATELLITE * 00000000 0 

FILTER STATE (EPOCH) —0.00809398 O 00226794-0.^0271521 0.00000625 O.OOOOOOIS 0.00000864 O OOOOOCMX1 0.00000000 0.00000000 
SMOOTHER STATE(EPOCH)*-0.01S3 163 1-0.00629638 0.07246483-0.00007387 0.00003446 0.00001248-0.00000000 0.00000000-0.00000000 
SMOOTHER SIGMAS(EPDCH) 0.00383015 0.01201316 O.OB267573 0.00007999 0.00003771 0.000002B3 0.00000000 0.00000000 C.OOOOOOOO 

NOMINAL TRAJECTORY (ECI)* 6211.3546-3075.3575 13992162 -1.80665 -033305 728347 
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SMOOTHER OUTPUT AT TIME O O 

GPAV GRAV 

O O 

O O 0.0 

0 . 00000000 - 0.00000000 
O OOOCOOOO 0.00000000 


NOMINAL TRAJECTORY (ECl)- 

-6413.0778 

2996.5868 

-0.6872 

0.45640 

0 97591 

-7 , 42596 

SMOOTHER 

TOTAL (ECI) 

-6413.0796 

2996.5888 

-0 €890 

0.45640 

0 97592 

-7.42596 

FILTER ESTIM. (HCL) 

0.0 

0.0 

0.0 

0.0 

0.0 

O O 

SMOOTHER 

ESTIM. (HCL) 

0.0023729 

Q.0007S36 

0.0020050 

0.0000030 

0.0000038 

-O 0000033 

SMOOTHER 

SIGMA (HCL) 

0.0040975 

0.0027577 

0.0125615 

0.0000120 

0.0000029 

0.0000043 


PARAMETER LABELS - X Y 2 XOOT TOOT 200T GTAV 

STATiON/SATELLlTE - O O O O O O O 

FILTER STATE (EFOCHl - 0.0 0.0 0.0 O O 0.0 O O 0.0 

SMOOTHER STATE (EPOCH)— 0.00 IS 1402 0, 00199352-0. 001T3A60-0.000001 1 1 0.00000421 0.00000331-0,00000000 
SMOOTHER SIGMAS(EPOCH) 0.0039C1tC 0.00372853 0.01237418 0.00001099 0.00000587 0.00000394 U.OOOOOOOO 


O 


o 
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FILTER POSITION SIGMA |flCL> 

• REPRESENTS 1C IN CCCI). H IN (HCU 

- REPRESENTS V IN |ECI1. C IN (IKL) 

♦ REPRESENTS 2 IN (ECi), fi. |N (HCL5 

-I , ^ 1 ^ 1 .., I 
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riLIER POSITION SIGMA (HCt) 

• REPRESENIS X IN (ECl). H IN (MCL| 

•> represents V IN (ECI). C IN CHCL} 

♦ represents 7 IN CECI>, L IN IHCU 
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FILTER POSITION SIGMA OlCL I 
• REPRESENTS X IN (ECl). H IM (IXIll 

- REPRESENfS V IN (ECl). C IN OCL) 

*■ REPRESENfS Z IN rECO L IN (ItCL) 


•173E-Ot ♦ 


o o 
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FILTER POSITION SIGMA (IKTU 

• REPRESENTS X IN CtCI). H IN CHCL) 

- REPRESENTS V IN lECI). C IN (HCLT 

♦ REPRESENTS 2 IN TECTl. L IN THCL) 


_I ♦ 1- 
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• -- MEAN- .260I1597E-01. SO- .45162 126E-01. RMS- .520584O2E-O1 
fUSED 332 OF 332 POINTS WITHIN SIGMA) 
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ORIGINAL' PAGE IS 
OF POOR QUALITY 
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MEAN- 
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SHOOlHf* POSITIOM SICMS CHCL» 

• REPRESENIS K IN lECII. H IN IHCL) 

- REPRESENIS V IN (ECO. C IN (HCt) 

♦ REPRESENIS Z IN (ECl). t |N luctj 


So 
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• — ,a9B74777E-02. SD» . 10701504E-02* RUS- . 3I72B2 14E-02 


lUSEO 332 OF 332 POINTS WITHIN 10. SIGMA) 
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_l 1 — - 
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SMOOTHER position SIGMAS (MCLl 

» REPRESENTS >C IN CEClK H IN <HCL> 

- REPRESENTS V IN (ECU, C IN (MCE) 

♦ REPRESENTS 2 IN (ECU. L IN (HCLl 

|, — 1 ^ * 1 ♦ 1 - 1 ♦ 1 ♦ 
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SMOOTHER POSITION SIGMAS IHCL) 

• REPRESENTS X IN (EC1)« H IN (HCL» 

- REPRESENTS Y IN UCIK C IN fHCl ) 

» REPRESENTS 2 IN (EC|), L IN (HCL) 
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6.0 JOB CONTROL LANGUAGE ORIGINAL PAGE IS 

OF POOR QUALITY 

The listing following this paragraph Is an example of the JCL required 
to run PREFER on the IBM 360/95, The core required to run PREFER (without 
overlays) using the loader Is approximately 370K. This will vary somewhat 
depending upon the program options (e.g. GPSPAC/TDRS measurements, plots, 
smoothed covariance file) and number of I/O buffers. The use of overlays 
would only reduce the storage by approximately 50K. When tracking from 
many GPS satellites Is processed, 1t may be necessary to Increase the 
region to BOOK. Printout from the job will define (approximately) the 
actual region used by PREFER. 


.// EXEC L0ADER,REi3I0N.i30-370K.PARri*-'SIZE*3700OO.EP«MAIN'' 

//SYSLIN DD 

/ / DD DSN*«<?<OB J , D I SP* < OLD » DELETE ) 

/ / DD DSN* YCBPC . i3SP ACE . OBJ . D I SP-SHR 

/ / DD DSN- YCBPG . CDC I 0 . OBJ , D I SP-SHR 

/ /GO . FT09P00 1 DD UN I T-2400-9 , VOL-SER-36493 , LABEL* < 3 r NL > . I N ) , 

/ / DCB- ( RECFM-VBS . LRECL*7204 . BLKS 1 2E-720S t BUFNO- 1 ) , D I SP- < OLD . KEEP ) 
//GO.FTlOFOOl DD UNIT-DISK, DSN-«<.!<MEAS7 D ISP- (NEW. DELETE) , 

// DCB- (RECFM-VBS, LRECL»72G*f, BLKSIZE-720S, BUFNO-1 ) , SPACE-(CYL, (2, 1 ) ) 
//GO.FTllFOOl DD DUMMY, DCB-BLKSIZE-100 

//GO. FT20F001 DD DSN-ORBIT. GTDS. SLP 19)50. DATA, DISP*SHR, DCB-BUFNO-1 
//G0.FT21F001 DD DSN-ORBIT. GTDS. T I MCOF. DATA, D I SP-SHR 
//GO. FT22F001 DD DSN-ORBIT. GTDS. SLPTOD. DATA, D I SP-SHR, DCB-BUFNO-1 
/ /GO . FT30FOO 1 DD UN I T-D I SK , DSN-«<2<U30 , D I S^» ( NEW, DELETE ) , 

// DCB- ( RECFM-VBS, LRECL-1 OS, BLKS I ZE-32 44, BUFNO- I ), SPACE- <TRK,( 10,2) ) 

//GO . FT40F00 1 DD DUMMY , DCB-BLKS I ZE* 100 

//G0.FT41F001 DD DUMMY, DCB-w J<SIZE*100 

//G0.FT42F001 DD DUMMY, DCB-BLKSIZE-100 

//G0.FT43F001 DD DUMMY, DCB-BLKSIZE-100 

//GO. FT50F001 DD UNIT-AFF-FT09F001 , LABEL»(4, NL, , IN) , VOL-SER-36493, 

/ / DCB- ( RECFM-VBS , LRECL-6-i.,S4 , BLKS I ZE-666S ) , D ISP- < OLD, KEEP ) 

/ /GO . M ASSTORE DD UN I T-D I SK , DSN-.V.?<U60 , D I SP- < NEW , DELETE ) , 

// SPACE-(CYL, (5,2) ) 

//GO. FT70F001 DD UNIT-DISK, DSN-.?<.1<0RB 1 , SPACE»(TRK, ( 10,2) ) , 

// DCB- ( RECFM-VBS, LRECL-X,BLKSIZE-7294, BUFNO- I ), DISP*( NEW, PASS ) 

//GO. FT71F00 1 DD UNIT-DISK, DSN-.V,S<TPH , SPACE- ( CYL, (2,2)), 

/ / DCB- (RECFM-VBS, LRECL-X • ’)^LKSI ZE-7294 > , DISP- ( NEW, DELETE ) 

//G0.DATA5 DD # 


Figure 7.1 Example of PREFER JCL for Processing Ground Tracking 
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Notice that Fortran units 9 and 50 use the same tape and tape drive. 
Unit 9 refers to file 1 (the observation file) and unit 50 refers to file 2 
(the ORBIT file). PREFER internally copies unit 9 to unit 10 (disk) so 
that the input files from the GTDS run may be placed on the same tape. 

Fortran unit MASSTORE is the scratch disk used by the filter and 
smoother to store the intermediate covariance matrices and other arrays. 

The DCB information is internal to PREFER. The space allocation of (5,2) 
cylinders should be adequate for most jobs. Unfortunately, it is very 
difficult to compute the storage required in a given run. It depends upon 
the number of states adjusted at each point in time. Since the number of 
states varies greatly during the run (as ground stations and GPS satellites 
change visibility to the satellite), the storage requirements can also vary 
greatly. In general, it should not be necessary to change the space 
allocation unless many GPS satellites are included. 

Fortran unit 30 is also scratch storage. However, the storage 
requirements are modest since only the filtered and smoothed state vectors 
are stored here. 

Units 70 and 71 are used for temporary storage of plotting variables 
and thus they cannot be ignored (e.g., "DUMMY") if plots are requested. 
Again the storage requirements are modest. Unit 71 is also used to output 
the smoothed covariance when this option is requested. The storage 
requirements for the smoothed covariance can be as large as 50 megabytes 
although 5 megabytes is a more typical number. 
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APPENPIX A 

PREFER Debug Printout 


FILTER (IBUG(D) 

NO - number of dynamic parameters 
TL - time of last entry to FILTER (sec) 

T - current time (epoch of mini-batch) (sec) 

X - state vector at epoch time of mini -batch 

"'*^**‘*x (TL to T) 

X - state vector (at time T) is printed after each measurement 
MEAS (IBU6(2)) 

TI - time of measurement (referenced to epoch of Host Trajectory) 

ISTA - station/ satellite number 

MTYPE - measurement type 

Y - measurement 

YD - measurement residual 

SIGM - input a of measurement 

PY - partial derivatives of measurement WRT state 

SNOISE (IBUG(2) 

T - current time (sec) 

OLDT - time of last entry to SNOISE (sec) 

NO - number of dynamic parameters 

X - state vector at time T (only dynamic states are printed) 

PHI - state transition matrix (OLDT to T) 

F - 6 X NP6 partial array (current cartesian 
elements WRT state vector) multiplied by DT 

A « F • QA where QA is spectral density matrix (state noise) 

Q * NO X NO state noise matrix in ECI coordinates 
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After transformation to epoch osculating elements, Q is printed 
again along with PHI 4 (the transition matrix from epoch osculating 
elements to current cartesian elements) 


DYNAM (IBUG(4)) 

TSTART - reference time (start of integration) 

TEND - end time of integration 
DT - TENO-TSTART 

STMAX - maximum allowed step size (with sign) for Taylor series 
integration 

HSTEP - current step size of interated integration 
HLFSTP »■ HSTEP/2 

TK - current time of iterated integration 
OLDTK - last TK 

PHI! - 6x6 transition matrix - cartesian (OLDTK) to cartesian (TK) 

PHI 4 - 6x6 transition matrix - epoch osculating to cartesian (TEND) 

PHI - NDxND output state transition' matrix 

XI - nominal cartesian elements at TEND 

XOUT - filter estimate of dynamic states at TEND 

CARD (IBU6(3) 

NUMCDS - number of adjusted parameters initialized 

LAB • parameter labels as read in 

SIG - a priori cr ' s as read in 

QN - state spectral density as read in 

STAVAR - station measurement variances as read in 

NORAG - number of drag parameters 
NGRAV - number of gravitational parameters 
NTHRST - number of thrust parameters 
NCLOCK - number of clock parameters 
HALT - number of altimeter parameters 
NPARM » NDRAG+3*(NTHRST+NGRAV) 
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NP6 » 6 + NPARM 

NO “ number of dynamic parameters (NP6+NCL0CK+NALT) 

NBIAS - 0 or 1 (1 indicates that biases are adjusted) 

NREFR - 0 or 1 (1 indicates that refraction parameters are adjusted) 

INTERG (IBU6(6)) 

NSAT - satellite number (1-24) 

TKl - current time in seconds from epoch of host tape 
TK - current time in seconds from epoch of GPS tape 
E - interpolated position and velocity at TK 

GETVEC (IBUG(6)) 

REQTIM - interpolation time 

XTO - interpolated satellite position and velocity at REQTIM 
EARTH (IBUG(7)) 

T - time from epoch 

STAPOS - station positions (Earth Centered Fixed) 

XSTA - station positions and velocities (Earth Centered Inertial) 

EVAL (IBUG(7)) 

CETOL - time tolerance for ephemeris computation 

NWSLP - Fortran unit number for SLP file 

REFDA2 modified Julian date corresponding to January 1950.0 

REFDA4 modified Julian date corresponding to initial conditions 

NCFOAY - number of days per curve fit 

TSEC - time in seconds from start of year to midpoint of this record time 
i nterval 

POELH - polynomial coefficients for delta H 

IDAY - beginning day of current record 

IN013 - polar motion switch (1 - compute polar motion) 
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USUN(3) - sun vector 

A - transformation matrix: selenocentric to selenographic 
ADOT - A 

B - transformation matrix: earth inertial (TOD) to earth fixed 
C - transformation matrix: mean of 1950.0 to true of date 
6HA - Greenwich hour angle 
XP - X polar motion angle 
YP - Y polar motion angle 

TZERO - number of A.1 seconds from 1950 to epoch 

T2ER01 - time in seconds from beginnning of ephemeris year to epoch 

OEPOCH - Julian date of epoch 

AZERO - difference between A.l and UTC time at epoch 
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APPENDIX B 

Guidelines for Setting Up GTDS Runs 

Since PREFER must use files created by GTDS, it is important that the 
GTOS run be set up correctly. The following guidelines are given. 

1) The GTDS run must create an Observation Savfi File (Unit 46) and 
an ORBIT file with partial derivatives (Unit 21) on magnetic tape. 
The input cards required to do this are: 


1 - 8 

9-11 

12 - 14 15 - 17 18 - 38 

39 - 59 

60 - 80 

SAVE 

1 



time 

OUTOPT 

0 

1 0 start time 

end time 

interval 


It is suggested that a time interval of 30 seconds be used. 

a • * * 

2) It is strongly recommended that the Observation Save File and the 
Orbit File be placed on the same tape (with the Observation file 
first) and tape drive (use UNIT=AFF). This will minimize the num- 
ber of tapes and drives required for both the GTDS and PREFER runs. 

3) Time regularized integration must not be used. 

4) Partial derivatives for the ORBIT file should be cartesian orbital 
elements with respect to epoch orbital elements. The epoch ele- 
ments may be cartesian, Keplerian or spherical, but cartesian are 
preferred. 

5) Since PREFER is intended to produce accurate ephemerides, it is 
important that the best models available be used in GTDS. In 
particular, the refraction corrections should be made in GTDS even 
though this increases the computer costs of a differential correc- 
tion run. This is done using an OBSCORR card. 
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APPENDIX C 
Subroutine GSPACE 

The subroutine GSPACE gives the user the capability to dynamically 
allocate core, that is, to allocate core for an array or a number of arrays 
during the execution of the program. This frees the user from hard^coding 
the dimension of arrays which will vary in size from run to run. 

GSPACE works by taking the amount of core (in bytes) requested by the 
user, incrementing it up to the next multiple of 8, and then issues a 
6ETHAIN instruction for that amount of core. GSPACE then calculates the 
offset, in numbers of 4 byte words, from some specific (hard) dimensioned 
variable to the first location of the space obtained. The user may then 
access this space by specifying the offset as an index to the (hard) 
dimension variable. 

General Format 

CALL GSPACE (N, NA, NB, SY, SZ) 


where 

N » the number of arrays for which core is being requested 
(INPUT TO GSPACE) 

NA » an integer array of dimension M; the 1st element of A is set 
equal to the amount of core (bytes) requested for the 1st 
array, the 2nd element is set equal to the amount of core 
requested for the second array, etc. (INPUT TO GSPACE) 

NB * an integer array of dimension N: the 1st element is set 
equal to the offset in 4 byte words from NA(1) to the 
storage area for the first array, the 2nd element is set 
equal to the offset in 4 byte words from NA(1) to the 
storage area for the second array, etc. (OUTPUT FROM GSPACE) 
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Y « the line to which control is transferred if the amount of 
core requested by the user is not available. (OUTPUT FROM 
GSPACE) 

Z ■ the line to which control is transferred if the call format 
is incorrect or the input variables contain invalid values 
(negative or real numbers) (OUTPUT FROM GSPACE) 

Example: 

niMENSION NA(2), NB(2) 

NA(1) - NPARMS * 4 
NA(2) » NOEV * 8 

CAU. GSPACE (2, NA, NB, « 100, A 200) 

To access the storage area for the 1st array, specify NA(NB(1) + I), I 
» 1, NPARMS. NPARMS 4-byte words can be stored in this space. To access 
the storage area for the 2nd array, specify NA(NB(2)+I), I = 1, N0EV*2. 
Either NDEV*2 4-byte words or NOEV 8-byte words can be stored in this 
area. Note that the offset to the 2nd storage area will be calculated from 
NA(1) and not NA(2). 

How GSPACE works: 

GSPACE first takes the input array NA values and checks to see if they 
are valid numbers. If they are, it then rounds them up to the next 
multiple of eight, even if they are already a multiple of eight. GSPACE 
then calls a routine which issues a separate GETMAIN instruction for each 
array requested. This will not necessarily result in the different array 
storage areas (within the same call to GSPACE) being contiguous (if the 
user has indeed requested multiple array storage areas). If the core is 
not available for any one of the GETMAIN requests, control is transferred 
to the line specified for Y. If all the core requested is available, 

GSPACE then enters a loop to calculate the offsets by the following 
algorithm: 
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ORIGINAL RMaEfis 
^FfRQGRCQlfAiTTY 


/storage address of first location 
NB(I) ■ array I storage area 


ys 


storage address' 
KA(1) 




When the loop Is complete, control returns to the calling program. Because 
the GETMAIN always gets core In the highest location (within the region) 
that is possible, subsequent GETMAIN's will likely get core closer to the 
load module (stored in the lowest region locations of core), and therefore 
the increments will decrease, i.e., NB(1)>NR(2)>NB(3). . . This will not 
necessarily be true if less core is requested by the subsequent GETMAIN's, 
i.e., NA(1)>NA(2)>NA(3), . . • 


REGION AVAILABLE 


NA(1) NA(2) 



/ 

/ 

/ 

1 


ARRAY 

2 

I/O 

BLOCKS 

LU 

QC. 

O 

O 

LU 

LU 

QC 

Li. 

ARRAY 

2 

Load Module 


NB(1) Offset 
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Other Entry Points 

CALL RELESE (N, NA, NB) - releases core allocated in a previous call 
to 6SPACE. N, NA and NB have the same definitions as for SSPACE. 

CALL REGION - computes IREG In common block C0RE$$. 

CALL GCORE - computes MCOR and ICORE 1n common block C0RE$S. 

Note; GCORE will allocate the entire region specified on 
the EXEC statement. Thus the user should Ignore 
the SIZE listed In the job system messages. 

Common Block/C0RE$$/ 


IREG - allocated region (In bytes) of current job step (from REGION 
parameter on EXEC statement). 

MCOR - maximum unused region (bytes) in any previous call to GCORE 
(including the current call). 


ICORE - unused region (bytes) at the current call to GCORE. 
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APPENDIX D 

Description of XOC Type" Disk 
Random Access I/O Package 


PREFER uses an assembly language subroutine which simulates "CDC type" 
random access disk I/O on an IBM 360. This subroutine has several 
advantages compared to standard IBM random access (or OAIO) 1n that the 
length and number of records can be variable. In fact, it Is possible to 
write over existing records with a d1 f ferent length as long as that length 
Is less than or equal to the old length. Although this routine is not 
suitable for use with pei'manent data sets, It 1s very useful for scratch 
storage within a program. This routine Is documented here since It is 
general purpose and may be useful for a variety of other applications. 

The subroutine has six entry points: OPENMS, ENDMS, WRITSS, WRITMS, 
REAOMS AND CLOSMS. A description of each of these functions is given below, 

1) CALL OPENMS (IX) - Opens mass storage (DONAME-MASSTORE) for write/ 
print /sequential access and saves the address of the record index for 
use by the routines WRITSS, REAOMS, and WRITMS, iX is an integer 
array which must be dimensioned as large as the number of records to 
be written. Note that the DCS information for MASSTORE is obtained 
internally in the routine. The user need only specify UNIT, DISP and 
SPACE. 

2) CALL WRITSS(IA,IR) - Writes a record to the next available space and 

stores the relative address in the index. If the entire record will 
not fit on the currant track, then the remainder of the record is 
segmented and written to subsequent tracks. lA is the data array to 
be written to disk and IR is the record number. Note : The first 

4-byte i nteger of lA must contain the length of the record in words 
(4-bytes). The example at the end demonstrates the use of WRITSS and 
READMS for writing and reading double precision words. 
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3) CALL ENOMS - closes the mass storage data set for sequential access, 
alters the macro form of the DCB to read/point and reopens the data 
set for update. 

4) GALL READMS(IA,IR) - reads a record from mass storage. If the record 
spans several tracks then it will be brought 1n segment by segment. 

The point macro Is used to position the disk at the beginning of the 
records. Each read attempts to read an entire track of data 
(Incorrect length Is ignored). The actual length of segments read in 
Is computed from the residual count In the CSW, and segments are read 
In until the total number of bytes read equals the record length. It 
Is assumed that the record length Is stored In the 1st fullword of the 
record* lA Is the array to which data will be retrieved from the disk 
and IR Is the record number. 

5) CALL WRITMS(IA,IR) - writes an updated record back to mass storage. 

If the record spans several tracks, then 1t will be rewritten in 
segments. The point macro Is used to specify the disk address of each 
segment. On each write, a length of one track is specified, and 
Incorrect length errors are intercepted by the SYNAO routine 'EX'. 

The actual length of a written segment Is determined from the residual 
count In the CSW. Segments are written until the total number of 
bytes written equals the record length. The definitions of lA and IR 
are the same as for WRITSS. Note that the length of the updated 
record must be less than or equal to the length of the old record. 

6} CALL CLOSMS - closes the mass storage dataset. 

ERROR CODES 

The routines can abend with 5 different error codes. Core dumps are 

produced by these abends. 
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di^lOfliVAL m 
OF Poor c- ^fy 

103 - error in WRITSS: too many bytes have been written. 

104 - error in READMS: too many bytes have been read. Record length 

was probably Incorrect. 

105 - error in WRITMS; too many bytes have been written. 

106 - error in WRITMS: I/O error other than incorrect length. 

999 - error in READMS: ENOMS was not called prior to first call to 

READMS. 



